Integral and Series Representations of Riemann's Zeta 
function, Dirichelet's Eta Function and a Medley of Related 

Results 

Michael S. Milgram * 

Consulting Physicist, Geometries Unlimited, Ltd. 
Box 1484, Deep River, Ont. Canada. KOJ IPO 

5: August 13, 2012 



Abstract 

Abstract: Contour integral representations for Riemann's Zeta function and Dirichelet's 
Eta (alternating Zeta) function are presented and investigated. These representations flow 
naturally from methods developed in the 1800's, but somehow they do not appear in the 
standard reference summaries, textbooks or literature. Using these representations as a basis, 
alternate derivations of known series and integral representations for the Zeta and Eta 
function are obtained on a unified basis that differs from the textbook approach, and results 
are developed that appear to be new. 

1 Introduction 

Riemann's Zeta function C,{s), and its sibling, Dirichelet's (alternating zeta) function •q{s) play an 
important role in physics, complex analysis and number theory, and have been studied extensively 
for several centuries. In the same vein, the general importance of a contour integral representation 
for any function has been knovifn for almost two centuries, so it is surprising that contour integral 
representations for both C,{s) and ry(s) exist that cannot be found in any of the modern handbooks 
(NIST, [29], Section 25.5(iii); Abramowitz and Stegun, [1], Chapter 23), textbooks (Apostol, [5], 
Chapter 12; Olver, [30], Chapter 8.2; Titchmarsh, [37], Chapter 4; Whittaker and Watson, [39], 
Section (13.13)), summaries (Edwards, [12], Ivic, [19], Chapter 4; Patterson, [31]; Srivastava and 
Choi, [34] and [35]), compendia (Bateman, [14], Section 1.12 and Chapter 17), tables (Gradshteyn 
and Ryzhik, [17], 9.512; Prudnikov et. al, [32], Appendix II.7) and websites ([13], [29], [40], [41], 
[42], [43]) that summarize what is known about these functions. One can only conclude that such 
representations, generally discussed in the literature of the late 1800's, have been long-buried and 
their significance overlooked by modern scholars. 



It is the purpose of this work to disinter these representations, revisit and explore some of the 
consequences. In particular, it is possible to obtain series and integral representations of C(s) and 
77(5) that unify and generalize well-known results in various limits and combinations and 
reproduce results that are only now being discovered and investigated by alternate means. 
Additionally, it is possible to explore the properties of C,{s) in the complex plane and on the 
critical line, and obtain results that appear to be new. Many of the results obtained are disparate 
and difficult to categorize in a unified manner, but share the common theme that they are all 
somehow obtained from a study of the revived integral representations. That is the unifying 
theme of this work. To maintain a semblance of brevity, most of the derivations are only sketched, 
with citations sufficient enough to allow the reader to reproduce any result for her/himself. 

* mike@geometrics-unlimited.com 
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2 Contour integral representations for ({s) and t]{s) 

Throughout, I shall use s = a + ip, where s G € (complex), a,pG^R (reals) and n and A'' are 
positive integers. 

The Riemann Zeta function ({s) and the Dirichelet alternating zeta function ry(s) are well-known 
and defined by (convergent) series representations: 



C(s) = ^fc-^ a>l (2.1) 
fe=i 

oo 

r?(s)^-^(-l)'=fc-« a>0 (2.2) 



fe=i 



with 



,7(s) = (l-2i-)C(5). (2.3) 

It is easily found by an elementary application of the residue theorem, that the following 
reproduces (2.1) 

C(s) = 77 residue I t=fc (p —— — - dt a > 1 (2.4) 
^ Mk sin(7ri) 

and the following reproduces (2.2) 



oo „ 

r]{s) = n'^residue\t=k f cot(7rt) e'''*dt ct > (2.5) 

fe=i 

where the contour of each of the integrals encloses the pole at t = k. 

Convert the result (2.4) into a contour integral enclosing the positive integers on the t-axis in a 
clockwise direction, and, provided that u > 1 so that contributions from infinity vanish, the 
contour may then be opened such that it stretches vertically in the complex t-plane, giving 

C(s) = - / coi{TTt)dt 0<c<l,ceai (2.6) 

^ J c—ioo 

1 r°° 

= --/ {c + it)~^ coi{-K{c + it))dt a > 1, (2.7) 

2 J-oo 

a result that cannot be found in any of the modern reference works cited (in particular [36]). It 
should be noted that this procedure is hardly novel - it is often employed in an educational 
context (e.g. in the special case s = 2) to demonstrate the utility of complex integration to 
evaluate special sums ([34] Section 4.1), ([15], Lemma 2.1), and forms the basis for an entire 
branch of physics ([27]). 

A similar application to (2.5), with some trivial trigonometric simplification, gives 



. .c+ioo ^ 



r](s) = - / — — - dt (2.8) 

' 2y,_,^ sin(7rf) ^ ' 

+ dt. (2.9) 



1 /-°° (cH 

2 sin(7r 



{c + it)) 

Eq. (2.9) reduces, with c = 1/2, to a Jensen result ([20]), presented in 1895, and to the best of my 

knowledge, reproduced, again, but only in the special case c = 1/2, only once in the modern 
reference literature ([41]). In particular, although Lindelof's 1905 work [22] arrives at a result 



2 



equivalent to that obtained here, attributing the technique to Cauchy (1827), it is performed in a 
general context, so that the specific results (2.7) and (2.9) never appear, particularly in that 
section of the work devoted to C{s). In modern notation, for a meromorphic function f{z), 
Lindelof writes ([22], Eq. /// (4)) in general 

oo If 

E f{'^) = -7r-f^^°^i^ ^) /(^) dz, (2.10) 



7y— — oo 



where the contour of integration encloses all the singularities of the integrand, but never 
specifically identifies f{z) = z~^. Much later, after the integral has been split into two by the 
invocation of identities for cot{'K z), he docs make this identification ([22], Eq. (6)), arriving 

(with 2.3) at (3.7 - see below) a now-well-known result, thereby bypassing the contour integral 
representation (2.7). In a similar vein, Olver ([30], page 290) reproduces (2.10), employing it to 
evaluate remainder terms, the Abcl-Plana formula and a Jensen result, but again, never writes 
(2.7) explicitly. Perhaps the omission of an explicit statement of (2.7) and (2.9) from Chapter IV 
of Lindelof ([22]) is why these forms seem to have vanished from the historical record, although 
citations to Jensen [20] and Lindelof [22] abound in modern summaries^. 

At this point, it is worthwhile to quote several related integral representations that define some of 
the relevant and important functions of complex analysis, starting with Riemann's original 1859 
contour integral representation of the zeta function ([33]) 

where the contour of integration encloses the negative t-axis, looping from t = — oo — zO to 
t = —oo + iO enclosing the point t=0. 

The equivalence of (2.11) and (2.1) is well-described in texts (e.g. [35]) and is obtained by 
reducing three components of the contour in (2.11). A different analysis is possible however (e.g. 
[12], Section 1.6) - open and translate the contour in (2.11) such that it lies vertically to the right 
of the origin in the complex t-plane^ with a <0 and evaluate the residues of each of the poles of 
the integrand lying a.t t = ± 2nni to find 

oo 

C(s) = — r(l - s)(27r«)^(exp(i7rs/2) - exp(-i7rs/2)) V k"''^. (2.12) 
This can be further reduced, giving the well-known functional equation for ({s) 

= (2,3) 

valid for all s by analytic continuation from a < 0. This demonstrates that (2.11) and (2.6) are 
fundamentally different representations of C{s), since further deformations of (2.11) do not lead to 
(2.7). 

Ephemerally, it has been claimed [41] that the following complex integral representation (missing 
a factor 27r in [41]) 



C(^) = T77r^/ dt a>0 (2.14) 



2{s-l)J_^ cosh2(7ri) 

can be obtained by the invocation of the Cauchy-Schlomilch transformation [4]. In fact, the 
representation (2.14) follows by the simple expedient of integrating (2.7) by parts in the special 
case c = 1/2 (see (3.16) below), as also noted in [41]. Another result 



^Recently, Srivastava and Choi [35], pages 169-172 have reproduced Lindelof's analysis, also failing to arrive 
explicitly at (2.7); although they do write a general form of (2.8) ([35], Section 2.3, Eq. (38)), it is studied only for 
the case c = 1/2. 

^and thereby vanishes 
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(i_2i-«)r(s + i)yo cosh2(t) 

obtained in [4] on the basis of the Cauchy-Schlomilch transformation does not follow from (2.14) 

as erroneously claimed in [41]. 

Based on (2.11), a number of other representations of ({s) are well-known. From the tables cited, 
the following are worth noting here: 
([14], Eq.l.l2(6)) 



([14], Eq.l.l2(4)) 

r(s) Jo exp(t) - 1 



1 roo 

C(s) = — - / — — -dt C7 > 1 , (2.17) 



C(s) = T^FTT^-TT / —TYT^'^t. a>l (2.18) 



and ([17], Eq. 3.527(1)) 

1 (2a)(''+i) f°° t 

4 r(s + 1) Jo ^M?(ai) 

By way of comparison, Laplace's representation ([17], Eq. 8.315(2)) 

,(c+it) 



1 _ 1 r 

T{s)~2^iJ_, 



:dt (2.19) 



{c + ity 

defines the inverse Gamma function. 

3 Simple reductions 
3.1 The case < c < 1 

The simplest reductions of (2.7) and (2.9) are obtained by uniting the two halves of their range 
corresponding to f < and < > 0. With the use of 

{-{c + i ty + {c-i ty) i = 2 (c2 + f2)(V2) sin(s arctan(t/c)) (3.1) 

and 

{c + ity + {c-i ty = 2 (c^ + i2)(«/2) cos(s arctan(i/c)) (3.2) 
one finds (with < c < 1 and < c < 1 respectively) 

(3.3) 

cos(s arctan(t/c)) cos(7r c) sin(7r c) — sinh(7r t) cosh(7r t) sin(s arctan(f/c)) 



(is) = -[ 
Jo 



(c2 +^2)^/2 (cosh2(7ri) - cos2(7r c)) 

a > 1 



-dt 



and 

^ .o.^/o^ ^2. o. ' ' ^-^dt, (3.4) 



-sin(sarctan(t/c)) cos(7rc) sinh(7rt) + cos(s arctan(f/c)) sin(7rc) cosh(7r<) 
(c2 + t^yli (cosh2(7rt) - cos2(7rc)) 



the latter valid for all s. In the case that c = |, (3.3) reduces to 

sinh(./)^u(.Hrc-|.a.(2f)) 
' Jo (l + 4i2)V2cosh(7rt) ' ^ ^ 
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equivalent to the well-known result"^ (which also appears as an exercise in many textbooks - e.g. 
[35], Eq.2.3(41)) 

C(a)-^-2^ /°° sin(^arctan(20)e-* 
^^'> - s-1 ^ io (1 + 4 t^y/' cosh(7r t) "'^ • 

Similarly, (3.4) reduces to the known ([22], IV.III(6)) result 

io (l + i2)V2cosh(7rt/2)"'*- ^"^-^^ 
Comparison of (3.5) and (3.6) immediately yields the known result ([17], Eq.3. 975(2)) 

^ 2(^-1) /"^ sin(sarctan(2i))e'^* , , 

^(^) = T37 + 2Vo (1 + 4.^)^2^4^) '^- ^'-'^ 

Notice that (3.6) and (3.7) are valid for all s. Simple substitution of c ^ in (3.4) gives 

??(s) = -sin(s7r/2) / — - — -dt a < (3.9) 



/o sinh(7r t) 
which, after taking (2.3) into account, leads to 



sin(s TT /2) 

(21-" - 1) Jq sinh(7ri) 



CG^) = i^^/ -TT^dt, a<0 (3.10) 



displaying the existence of the trivial zeros of ({s) at s = —2N, since the integral in (3.10) is 
clearly finite at each of these values of s. The result (3.10) is equivalent to the known result (2.16) 
after taking (2.13) into account. Other possibilities abound. For example, set c = 1/3 in (3.3) to 
obtain 

_ os-i (2 siii(s arctan(.)) sinh(2 tt t /3) - V3 cos(.s arctan(.))) 
' Jo (l+t2)«/2(2cosh(2^</3) + l) ■ ^' ' 

Further, a novel representation can be found by noting that (2.9) is invariant under the 
transformation of variables t ^ t + p together with the choice c = a along with the requirement 
that < o" < 1. This effectively allows one to choose; c = ,s, yielding the following complex 
representation for ri{s), as well as the the corresponding representation for ({s) because of (2.3), 
valid in the critical strip < u < 1: 



' r <• + "'"• (3.13) 



2 y_ooSin(7r (s + it)) 

Because (3.1) and (3.2) are valid for c G £, it follows that (3.4) with c = s is also a valid 
representation for r]{s), as well as the the corresponding representation for ({s) with recourse to 
(2.3), again only valid in the critical strip < < 1. Another interesting variation is the choice 
c = 1/p in (3.4), valid for p > 1, again applicable to the corresponding representation for ({s) 
because of (2.3). Such a representation may find application in the analysis of ^(0" + ip) in the 
asymptotic limit p ^ 00. Finally, the choice c = a, 0<a<l leads to an interesting variation 
that is discussed in Section 9. 



^attributed to Jensen by Lindelof, and a special case of Hermite's result for the Hurwitz zeta function ([39], 
Section 13.2). 
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3.2 Simple Recursion via Partial Differentiation 

As a function of c, ({s) and r]{s) are constant, as are their integral representations over a range 
corresponding to < N — c < 1; however the integrals change discontinuously at each new value 
of c that corresponds to a new value of TV. Thus, the operator ^ acting on (2.7) and (2.9) 
vanishes, except at non-negative integers c = N, when it becomes indefinite. Therefore, it is of 
interest to investigate this operation in the neighbourhood of c = N/2, as well as the limit N. 
Applied to (2.9)^, the requirement that 

dr]{s)/dc = (3.14) 

and identification of one of the resulting terms, yields 

2s J _^ sm (7r(c + tt)) 
a new result with improved convergence at infinity. Similarly^ , operating with 4- on (2.7) yields 



a.+i) = ^r .y/t' (3.16) 

2s y_ooSm^(7r {c + it}) 

Using (3.1) and (3.2) to reduce the integration range of (3.15), gives a fairly lengthy result for 
general values of c. If c = 1/2, (3.15) becomes 

2-. ^-(l + .^)-V^sinh(.t/2)sin(.arctan(t)) 
' s Jo cosh2(7ri/2) 
which, by letting s ^ s — 1 can be rewritten as 

,(. + l) = ^ii^M£)+?!:iZ[ rcos(sarctan(t))(l + .^)-V^sinh(.V2)^^^^ 
s s Jo cosh (7rt/2) 

and, with regard to (2.3) 



2(s-l)(1^2-^ 
^ s(l-2-«) 

2^-1 TT f°° cos(s arctan(t)) (1 + i2)-V2 sinh(7r ^2) 



(3.19) 



7 

Jo 



tdt . 



s {1-2-') Jo cosh^(7rt/2) 
Similarly^ , in the case that c = 1/2, (3.16) reduces to the well-known result ([22], Eq. 4,//7(5)'^) 

Cis + 1) = ^ r ^os{s^rctan(f^ 

' s Jo il + f^y/^cosh\nt/2)'^^- ^^-^"^ 

If c in (3.15)^ , one finds, a finite result^. 



, S TT 



7rsin( — ) .oo^l-. w ^^ 

V{s)= ^ ; / —r-j-^dt, a<0 3.21 

s-1 Jo smh (tt*) 



^and equivalent to integration by parts 

■'It should be noted that Srivastava, Adamchik and Glasscr ([36]) study an incomplete form of (3.16) in order to 
obtain representations for sums both involving C('^)> s-nd CC") by itself, witliout commenting on the significance of 
the contour integral with infinite bounds (3.16). 

®(3.18) and (3.19) are the first two entries in a family of recursive relationships that will be discussed elsewhere 
(in preparation). 

'^attributed to Jensen [20] 

^For the case c = in (3.16) see the following section 
^also equivalent to integrating (3.9) by parts 
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and, when c = 1, one finds the unusual forms 

iirr r{l + rie' ^)(-*) cosh(7r re'^)e'^ 



r TT f 

s J I sinh''(7rri) 



(3.22) 



sinh^(7rr e'^) 
r TT (1 + f2)(-s/2) (,og(-g arctan(r t)) cosh(7rrt) 



d6» 



or, 



r]{s + 1) = 1 + —- / ^ ^ '- de 

r TT /■°° (1 + t2)(-s/2) (,Qg(g arctan(r <)) cosh(7r r t) 



(3.23) 



sinh (-Trri) 



■dt 



valid for any^" < r < 1 and all s, after deformation of the contour of integration in (3.15) either 
to the left or right of the singularity at t = 1, respectively. 

4 Special cases 
4.1 The case c=0 

With respect to (3.16), setting c = together with some simplification gives 



C(S) = -TT^ 



,sin(7rs/2) t^-' , „ 



s-1 Jo sinh^(t) 

reproducing (2.18), with a= 1, after taking (2.13) into account. This same case (c = 0) naturally 
suggests that the integral in (2.6) be broken into two parts 

C(s) = Co(s) + Ci(s) (4.2) 

where 

Co(s) = ^ (f t~'coi(nt)dt , (4.3) 
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(^i(s) = - / cot{-Kt) dt + - cot{-Kt) dt 



(4.4) 



such that the contour in (4.3) bypasses the origin to the right. It is easy to show that the sum of 
the two integrals in (4.4) reduces to 

/oo 
coth.{-Kt)dt a>l (4.5) 

sin(7rs/2) . , f°° e^-""'^ , 

= — — + sin(7r s/2) / —— — - dt , (4.6) 
s - 1 ^ ' ' Ji sinh(7r t) ' ^ ' 

the latter valid for all s. One, of several commonly used possibilities of reducing (4.3) is to apply 
([18] (54.3.2)) 

2 °° 

t cot(7rt) = ^t^'' C(2fc) \t\ < 1 (4.7) 

giving 



^"particularly r = l/7rorr = l/2 



7 



fc=0 

convergent for all s. Such sums appear frequently in the literature (e.g. [8], [34] Section (3.2)), 
and have been studied extensively when s = —n ([34], [7]). In the case s = 2n, (4.8) reduces to 
Co(s) = C{s) and (4.6) vanishes. 

Along with (4.1) the integral in (3.16) can be similarly split, using the convergent series 
representation 

^ oo 

= (2^-1) k) v' \v\ < 1 (4.9) 

sin (ttv) ^ 

where the Bernoulli numbers B2k appearing in the original ([18], Eq.(50.6.10)) are related to 
C(2fc) by 



Thus the representation 



1 in r t^-" , 7rsin7r.s/2) t^-* , 

C(s) = — -/ — ^ dt i '—^ 5 dt (4.11) 

' 2 s-1 J-i sin2 TTt s-1 Ji sinh'(7rt) ^ ' 



becomes (formally) 



^ s-1 [t^ fr'o 2A;-s A sinh2(7ri) y ^ ' 

after straightforward integration of (4.9) with a < 0. In (4.12), the integral converges, and the 
formal sum conditionally converges or diverges according to several tests. Howver, further analysis 
is still possible. Following the method of [34], Section (3.3) applied to (4.8), together with [34], 
Eq. 3.4(15), we have the following general result for \v\ < 1 

fe=o -"J fc=l 

1 1 r 

= — + -v'' t-" {^{l+t)-^il-t))dt a<2 (4.14) 

2 s 2 Jq 

so that (formally), after premultiplying (4.14) with 1/v, then operating with v^-^, 



E C(2fe)(2 fc 1)^^^ =-^ + h' f^'-' (*'( 1 + *) + 1 - t)) dt (4.15) 



2k-s 2s 2 

k=0 



= -^ + lv' ft^-' (tt' (1 + coe{nt)) -hdt. a<2 (4.16) 

Notice that the right-hand sides of (4.14) provides an analytic continuation of the left-hand side 
for \v\ > 1 and the right-hand side of (4.16) gives a regularization of the left-hand side. Putting 

together (4.12) and (4.16), incorporating the first term of (4.16) into the second term of (4.12), 
setting V = i and constraining cr > eventually yields a new integral representation^^ 

-k' 'sm(— ) .oo ^ ^ 



icompaxe with [38], Eq. (2.8). 
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valid for < (7 < 2, a non-overlapping analytic continuation of (4.1) that spans the critical strip. 
Additionally, using (2.13) followed by replacing s ^ 1 — s gives a second representation 



valid for — 1 < ct < 1, a non-overlapping analytic continuation of (2.18). Over the critical strip 
< (7 < 1 representations (4.17) and (4.18) are simultaneously true and augment the rather short 
list of other representations that share this property ([6], as cited in [14]^^). See Section 9 for 
further application of this result. 

4.2 The case c > 1 

Interesting results can be obtained by translating the contour in (2.6) and (2.9) and adding the 
residues so-omitted. From (2.1) and (2.2) define the partial sums 

N 



Civ(s)^^A;-« (4.19) 
fe=i 

N 

w(s) = -l](-l)'^fc-^ (4.20) 



fe=i 

let cn,N = 1, 2, 3, ... with 

N<CN<N + 1 (4.21) 

and adjust^^ all the above results accordingly by redefining c—^Cn. This gives, for example, the 
generalization of (2.6) 

C(s) = Cjv(s) - i^Sn,c^ - - / t-' cotiwt) dt ,a>l (4.22) 

^ JcN-ioo 

and (2.9) 

nis) = V.is) + i^^.,.. + - y^^_^^ dt (4.23) 

so (3.3) and (3.4) with cn become integral representations of the remainder of the partial 
sums C(s) — C,n{s) and r}{s) — r]N{s), the former corresponding to a special case of the Hurwitz 
zeta function ({s, N). 

In the case that cn = N = 1, one finds 



and 



^, , 1 /■°°cosh(7ri)sin(sarctan(i)) , ^ 

_ s + l ^J-e--^^^^^^ ^^^^^ 



2(s-l) Jo (l-Fi2)V2sinh(7ri) 



, , 1 /'°° sin(sarctan(t)) ,^ ,^ 

with the latter two results valid for all s. Neither of (4.24) nor (4.26) appear in references; (4.25) 
corresponds to a limiting case of Hermite's representation of the Hurwitz zcta hmction ([14], Eq. 



l^unnumbered equations immediately following 1.12(1), also reproduced in [35], Eqs. 2.3(44)-(46). 

^^If cjv = ^ an extra half-residue is removed and the contour deformed to pass to the right of the point t = N. 
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1.10(7)). In analogy to the developments leading to (3.8), equations (4.24), (4.26) and (2.3), 
together with the following modified form of ([17], Eq.(9.51.5)) 

C(,)- '^"'^ I /■- sin(.arctan(^))(l + e-*) 

(2.-l)(s-l) + 2«-iyo (l + i2)V2sinh(7rt) ^ ' 

can be used to obtain the possibly new result 

s-3 /"^ sin(sarctan(i))e''* , , 

For c = 1/2 + corresponding to general partial sums, two possibilities are apparent. Based on 
(3.1), (3.2), (3.3) and (3.4)1^ we find 



^11 r 



sinh(7r i (1 + 2 A^)) (1 + t2)-s/2 sin(s arctan(f)) , 

at (T > 1 

cosh2(7rt (1/2 + 7V)) 

(4.29) 

(l/2 + 7V)i-« /■'^(l + i2)-V2sin(sarctan(t)) 



s-1 

and 



which together extend (3.5) and (3.7) respectively. Application of (3.15) and (3.16) with (3.1) and 
(3.2), also with c = 1/2 + A'', gives the remainders 

as + -J1I1±^ r (l + ^-);^-cos(^arctan(0) 

fe^ ^ -^0 cosh2(7rt(l/2 + 7V)) ^ ^ 



k=l 

and 



''(^ + i) + Eyr 
fc=i 



(4.33) 



TT (-1)^ (1/2 + Ny-' f°° sinh(7r t (1/2 + TV)) (1 + i2)-«/2 sin(s arctan(t)) 



/ 

Jo 



cosh2(7rt(l/2 + AT)) 

valid for all s. For the case c = N applying (4.22) and (4.23) gives the remainders 

C(.) - f 1 = - J- + N^- rcosh(.Art)(l + .^)-V^sin(.arctan(t)) 

2Ar* Jo sinh(7rAri) ^ ^ 

= -^ + ^+2Ar-^ r (l+^^)-/^sin(.arctan(^)) 

2iV- + s-l+ io e2^-*-l ^ 

and 

^ A;^ 2Ar« ^ ^ Jo sinh(7rAri) ^ ^ 

which extend (4.24) and (4.26), the latter valid for all s. It should be noted that (4.35), valid for 
s G € generalizes well-established results (e.g. [17], Eq. 3.411(12)) that are only valid in the 



*for comparison, see also [30], Section 8.3, and [12], Section (6.4) and Chapter 7. 
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limiting case s = n^^. Omitting the integral term from (4.30) gives an upper bound for the 
remainder; similarly, omitting the integral term in (4.35) gives a lower bound for the remainder, 
provided that ,s G in both cases. In addition, many of the above, in particular (4.30) and (4.35) 
appear to be new and should be compared to classical results such as [19], Theorem 1.8, page 23 
and pages 99-100. 



5 Series Representations 
5.1 By Parts 

By simple manipulation of the integral representations given, it is possible to obtain new series 
representations. Starting from the integral representation (3.10), integrate by parts, giving 

({s)= sm(7rs/2)7r^ 



(2(i--)-l)(-s/2+l) 
^^-iF.(-| + l; ^,2-|; ^)cosh(7rt) 



dt cr < 



sinh^(7r t) 

where the hypergeometric function arises because of [23], Eq. 6. 2. 11(6). Write the hypergeometric 
function as a (uniformly convergent) series, thereby permitting the sum and integral to be 
interchanged, apply a second integration by parts, and find 



( \ °° 2k foc ,l-s+2k 

JT^) FpIT2)i l^-^' 

Identify the integral in (5.2) using (2.18) giving 

, _ 2^ w'-^ sin(7r s/2) ^ 2-^>' ({-s + 1 + 2 k)r{-s + 2 + 2k) 
(,{s) - A. r(2fc + 2) ^^""^^ 

valid for all values of a by the ratio test and the principle of analytic continuation. This result, 
which appears to be new, can be rewritten by extracting the fc = term, and applying (2.13) to 
yield 

, _ 2^7r^-isin(7r5/2) ^ 2-^1^ Cjs + 1 + 2 k)Ti-s + 2 + 2 k) 

21--2 + S ^ r(2fc + 2) ' ^ ' 

fc=i ^ ' 

a representation, valid for all values of s, that is reminiscent of the following similar sum recently 
obtained by Tyagi and Holm ([38] Eq.(3.5)) 



7r°-- (2-2«-2'=)C(-s-|-l + 2fc)r(-s-|-l + 2fc) 
Cis) = T3^sm(. s/2) ^ i r(2fc + 2) ' ^>'' 

(5.5) 



r(2fc + 2) 



Eq. (5.4) can also be rewritten in the form 
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1 f. 2-^^C(^ + 2fc)r(l + 5 + 2fc) 

^~r(s)(2--l-s)^ r(2A; + 2) ^^'^^ 



-'^Although [17], Eq. 3.411(21) gives (with a missing minus sign) an integral representation for finite sums such 
as those appearing in (4.35), specified to be valid only for integer exponents (fe"), that particular result in [17] also 
appears to apply after the generalization n — > s. 

^^by Gauss' test 
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by replacing s — 1 — s and applying (2.13). 

With reference to (2.7), express the numerator cosine factor as a difference of exponentials each of 
which in turn is written as a convergent power series in tt (c + it), recognize that the resulting 
series can be interchanged with the integral and each term identified using (2.9) to find 

^ /^V^ (-^V4)'^??(.s--2fc)\ 

A more interesting variation of (5.7) can be obtained by applying (2.3) and explicitly identifying 
the fc = term of the sum giving 

a convergent representation if ct > 1 (by Gauss' test). An immediate consequence is the 
identification 

C(2) = -y C(0). (5.9) 

because only the k = 1 term in the sum contributes when s = 2. Additionally, when s = 2n, (5.8) 

terminates at n terms, becoming a recursion for C(2't.) in terms of C(27^ — 2fc), k = 1, . . .n - see 
section 8 for further discussion. If the functional equation (2.13) is used to extend its region of 
convergence, (5.8) becomes 

.y^ C(l-s + 2fc) (2^-^^-1-1) 
r(2A: + l)r(s-2fc) 

^(-)=2 "\2-^-l)cos(../2) « 

which augments similar series appearing in [34], Sections 4.2 and 4.3, and corresponds to a 
combination of entries in [18] (Eqs. (54.6.3) and (54.6.4)). 

5.2 By splitting 
5.2.1 The lower range 

After applying (2.13), the integral in (4.1) is frequently (e.g. [35]) and conveniently split in two, 
defining 

p UJ J.S 

C-{s)= -TjjrJt (5-11) 

^0 smh (t) 

and 



with 



Tis + l) 
2 

for arbitrary w. Similar to the derivation of (4.12), after application of (2.13) we have 



C(s)4l^=C-W + C+W (5.13) 



k=0 

which can be rewritten using (4.10) as 



^""h - 1/ 2) r(fc + 1) (2 fc - 1 + s) 



(5.14) 
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~ {-'^) r(t + i)c(2fc) 



U) < TV 



(5.15) 



Substitute (2.1) in (5.15), interchange the order of summation, and identify the resulting series to 
obtain 

= ^i^H^^ ^ 2 + 2^2 + 2^ -;^)-^-th(a;) + ^ , (5.16) 

a convergent representation that is vahd for all s. 

Other representations are easily obtained by applying any of the standard hypergeometric linear 
transforms as in (6.4) - see below. One such leads to 



^2 2^^r(fc+i + |)^^V2mVa;2 + l 



f+'- w^coth(w) + 



SUJ 



s-l 



1 



(5.17) 



a rapidly converging representation for small values of w, valid for all s. The inner sum of (5.17) 
can be written in several ways. By expanding the denominator as a series in l/{t^rin?) with 
r < 2p — 1 and interchanging the two series we find the general form 



E 

m=0 



T{p + m) (-l)™f-2m-2p^(-2^ + 2j9-r) 
r(p) r(m + 1) 



(5.18) 



where the right hand side continues the left if \t\ > 1 and r — 2p> —1. Alternatively, [18], Eq. 
(6.1.32) identifies 



E 

m=0 



1 1 TT -uf^y\ 

(m2 t^ + y^) " V ^ 27^^^°^ 



(5.19) 



so for p = fc + 1 and r = in (5.18), each term of the inner series (5.17) can be obtained 
analytically by taking the fc*'' derivative of the right-hand side with respect to the variable at 
y = 1, leading to a rational polynomial in coth(a;) and sinh(a;). The first few terms are given 
below for two choices of w. 



The first few values of «;(0, A; + 1, tt^/w^) in (5.18) 



analytic 



0J = 1 



u = n 



-1/2+ wcoth(a;)/2 



1^1 cv^ 
2 4 sinh^w 



H — w coth(a;) 
4 



3a; j.u3 1 iv2 , cjcothcjf — 3+2li; ) 

coth w — i coth uj H , y ^ 

I 8+3 ui^ 

~^ If) sinh^c^' 



0.156518 1.076674 
0.927424 X 10"^ 0.306837 

0.795491 X 10"'^ 0.1.34296 



Application of (5.18) to (5.17) leads to a less rapidly converging representation that can be 
expressed in terms of more familiar functions 
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. , ,^ (5-20) 



^2 2'^^^^ r(fc+| + |)r(m + l) 

— w^cothfw) -I . 

s — 1 

A second transform of the intermediate hypergeometric equation yields another useful result 
2sL0^-' ^ ^=+1)^ u.'^ „^(_l)™^(m + A;+i + |)C(2m + 2A; + 2) 

^-^'^-r(i/2 + 5/2)^/^^ ^/^^ r(m + i)r(fc + i)(2fc + i + s) 

„ , ,s-l 

— a;*coth(a;) 



s- 1 

(5.21) 

Both (5.20) and (5.21 ) can be rc-ordcrcd along a diagonal of the double sum-'^^, allowing one of the 
resulting (terminating) series to be evaluated in closed form. In either case, both eventually give 



oo 



2 



C-(«) = -2sa;-i^(-J^) ^(^^-^^coth(a;), (5.22) 

reducing to (4.8) in the case cj = tt taking (2.13) into consideration. Furthermore, by applying 
(2.1) and interchanging the resulting series, a further (more slowly converging) representation is 
found 

s — 1 ^ 1 -1 2 

m=l 

which differs from (5.16)^^. Series of the form (5.22) have been studied extensively ([34], [35]) 
when s = n for special values of w. The representations given above will find application in 
Section 8. 



5.2.2 The upper range 

Rewrite (5.12) in exponential form and expand the factor (1 — e~^*)~^ in a convergent power 
series, interchange the sum and integral and recognize that the resulting integral is a generalized 
exponential integral ^s{z) usually defined in terms of incomplete Gamma functions ([25], [29], 
Section 8.19 ) as follows 

Y,,{z) ^ z-'-^T {\- s,z) (5.24) 
giving a rapidly converging series representation 

oo 

C+(s) =4w"+^ ^/cE_^(2wfc). (5.25) 

k=l 

A useful variant of (5.25) is obtained by manipulating the well-known recurrence relation ([29], 
Eq. 8.19.12) for E^(0) to obtain 

E_,(z) = ^ (e-^(s + z) + s{s - 1) E2_«(z)) . (5.26) 

Apply (5.26) repeatedly A^" times to (5.25) and find 

"Er=o E^=o k) = Er=o E^=o k - m) 

^*and could have also been obtained by recourse to contiguity relations. 
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c,(.) = iog(e- - 1) + 2..- ^ + r(. + i)g 

2cj-''-^-^r(.s + l) ^ Ejv+2-.(2fca;) 
r(s - TV - 1) 2^+1^ fc^+i 

(5.27) 

in terms of logarithms and polylogarithms {Lik(z)), a result that diverges numerically as A'' 
increases, but which also, usefully, terminates if s = n, N > n > 0. This is discussed in more 
detail in Section 8. 

Another useful representation emerges from (5.25) by writing the fimction Es(..) in terms of 
confluent hypergeometric functions. From [25], Eq.(2.6b) and [21], Eq.(9.10.3) we find 

(5-28) 

C+(s) = 4a;«+i^fcE_,(2A;w) = 

fc=i 

A S-\- 1 ^ 

2i-« r(s + 1) C(s) - V fc 6-2'=- iFi(l; 2 + s;2ku), 

^ + ^ fe=i 

where the series on the right hand side converges absolutely if cr > 1 independent of u>. 
Comparison of (5.25), (5.13) and (5.28) identifies 

4 S+ 1 ^ 

C-(s) = — V&e-2'='"iFi(l; 2 + s; 2koj) (5.29) 

s + 1 

k=l 

which can be interpreted as a transformation of sums, or an analytic continuation, for the various 
guises of Q-{s) given in Section(5.2.1), if both, or either, of the two sides converge for some choice 
of the variables s and/or cj, rcspcictively. The fact that a term containing C,{s) naturally appears 
in the expression (5.28) for C+(s) suggests that there will be a severe cancellation of digits 
between the two terms C-(s) and (^+{s) if (5.13) is used to attempt a numerical evaluation of C,{s), 
as has been reported elsewhere ([11]). Essentially, (5.13) is numerically useless (with common 
choices of arithmetic precision) for any choices of p > 10, even if w is carefully chosen in an 
attempt to balance the cancellation of digits. However, with 20 digits of precision, using 20 terms 
in the series, it is possible to find the first zero of C(l/2 + ip) correct to 10 digits, deteriorating 
rapidly thereafter as p increases. For any choice of oj , the convergence properties of the two sums 
representing C,-{s) and C+(s) anti-correlate, as can be seen from the Table in the previous Section. 
This can also be seen by rewriting (5.28) to yield the combined representation 

2i-r(s + l)C(s)=4u;«+i (^A:(E_,(2fcu;) + i^^iii±li^^±^i^^)j . (5.30) 

In (5.30), the first inner term (E_s(2fca;)) vanishes like exp(— 2 fca;)/(2 fccj) whereas the second 
term (ii^i/(s + 1)) vanishes like l/fc(^+*) — exp(— 2 /c a;)/(2 few) so for large values of k cancellation 
of digits is bound to occur 

It should be noted that (5.28) is representative of a family of related transformations. Although 
setting s = in (5.28) does not lead to a new result, operating first with ^ followed by setting 
s = leads to a transformation of a related sum 

(5.31) 



a;arctanh(^^) ^a;arctanh(— |^) i ^ i ^ 
^Ei(2fcw) = r7=4^+9log(-) + ?77 + ^ 



k=l ^ k=2 



^^In fact, the coefficients of both asymptotic series prefaced by the exponential terms cancel exa<;tly term by term. 
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where the first term corresponding to fc = 1 has been isolated, since it contains the singularity 

_d_ 

ds 



that occurs due to divergence of the left-hand sum if lj = iw. Similarly, operating with 4- on 



(5.28) at s = 1, leads to the identification 

oo 

E + ^) = -1 - l°g(2) + log(e2 - 1) (5.32) 
fc=i ^ 

(which can alternatively be obtained by writing the infinite sum as a product). A second 
application of ^ on (5.28), again at s = 1 leads to the following transformation of sums 



fe=i 



(5.33) 



1^, ,7r2p + 1 , ^ 1 V-2 + 2e% ^Ei(2A;) 

k=l k=l 

This transformation can be further reduced by applying an integral representation given in [29] , 
Section 6.2, leading to 

Correcting a misprinted expression given in ([3], Eq. (4.8), [34] and [35], Eq. 3.3(43)) which 
should read 

^ ^ C(2 k) = log(7r csc(7r V^))—, (5.35) 
followed by a change of variables, eventually yields the representation 

f: ^'""-,'r-""' = -^/i(>OE(ie' - ie-') - M*)).*. (5.36) 

Substituting (5.32), (5.34)) and (5.36) into (5.33) followed by some simplification, finally gives 

.OS ( '^-'-y-"l ^ = + ^ + .(1) - i .»S-^(2) + 2 . (5.37) 



/ 

Jo 



2v J V 12 2 ' ' 2 



a result that appears to be new. 

5.3 Another Interpretation 

These expressions for the split representation can be interpreted differently. For the moment, let 
a < —3, and insert the contour integral representation [25], Eq. (2.6a) for E_s(z) into the 
converging series representation (5.25), noting that contributions from infinity vanish, giving 

a(.) = -4c.-«f::^ r"°"'-"'""''d. (5.38) 

s+i + t 

where the contour is defined by the real parameter c < —2. With this caveat, the scries in (5.38) 
converges, allowing the integration and summation to be interchanged, leading to the identification 

2 TT I Jc-i oo S+l+t 

The contour may now be deformed to enclose the singularity at t = — s — 1, translated to the right 
to pick up residues of r(— f) at t = n, n = 0, 1 . . . and the residue of C(~l — t) at t = —2. After 
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evaluation, each of the residues so obtained cancels one of the terms in the series expression (5.15) 
for (s) leaving (after a change of integration variables) 

r(. + i)C(.)_ 1 / r(. + i)CW 



the standard Cauchy representation of the left-hand side, valid for all s, where the contour 
encloses the simple pole at t = s. 



6 Integration by parts 

Integration by parts applied to previous results yields other representations. 
Prom (2.17), we obtain^" 

1 f°° 

'^('^ = - Y{s-l) I i'"'log(l-exp(-i))rft a>l. (6.1) 



and from (4.24) we find 

1 s /•°° log(sinh(7rt)) cos((s + 1) arctan(^)) 
^^'^=2-nJo (l + ^^)a+^)/^ ^^-'^ 

as well as the following intriguing result^^ 



1 1 V^^"^(^)r(^-1) /■o° tsiu(.sairUm(^)) 
C(sj =0 + TsZ + ^ / ■ .2/ 21*1(1/2, s/2; 3/2; ) dt 

2 2 p(_) 7o smh {-Kt) 



(6.3) 



/■°°fcosh(7rt)cos(sarctan(t)) ^ ,^ . „. „, , 

-W ■ u ^.2^ 2F1 1/2, a/2; 3/2; -t^ . a>0 

Jo smh(7rt) (1 +t^) 

This result has the interesting property that it is the only integral representation of which I am 
aware in which the independent variable s does not appear as an exponent. Many variations of 
(6.3) can be obtained through the use of any of the transformations of the hypergeometric 
function, for example 

2Fi(l/2, s/2; 3/2; -t^) = (1 + t')-^/^ 2Fi(l, s/2; 3/2; f/il + f)), (6.4) 
which choice happens to restore the variable s to its usual role as an exponent. 



7 Conversion to Integral Representations 

Ftom (5.4), identify the product C(...)r(...) in the form of an integral representation from one of 
the primary definitions of ({s) ([17], Eq. (9.511)) and interchange the scries and sum. The sums 
can now be explicitly evaluated, leaving an integral representation valid over an important range 
of s 



. / 2ssinh(V2) , , , A 

4^ sin(. s/ 2) r ' V 1 - ' + ) 

2 - 21+" + 5 2" Jo e* - 1 ^ ' 

Each of the terms in (7.1) correspond to known integrals ([17], Eqs. (3.552.1) and (8.312.2)), and 
if the identification is taken to completion, (7.1) reduces to the identity ({s) = ({s) with the help 



^"known to Maple, but only in the form of a Mellin transform 
discovered by the Maple computer code 
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of (2.13). Other useful representations are immediately available by evaluating the integrals in 
(7.1) corresponding to pairs of terms in the large parenthesis, giving 



7r^2«-2 /^sin(7rs)2^ /■°°* '(-l + ^o^M^)) l\ 
^(-)= (l-2^)cos(.a/2) |-^yo e^-1 ''^W)) ^'''^ 

and 

({s) = tt'-'^ sin(7r s/2) ( f u""-^ (w/sinh(t;) - 1) e"" dv - — — ) . a<2 . 

^ ^ ^ ' \Jo \ ! I sm(7rs)r(l + s); 

(7.3) 

8 The case s = n 

Several of the results obtained in previous sections reduce to helpful representations when s = n. 

8.1 Finite series representations 

From (5.4) we find, for even values of the argument, 

" (-l)fc7r^fc(2n-2fc^l)C(2n-2fc) 
ih r(2/c + 2) 

'^('-) = -2- F^^TTTT^l ' 

from (5.8) we find 

1 ^ (-7r^)^(l-2i-^"+^^)C(2n-2fc) 
C(2«)- 2(2-2" _ 1)1. r(2fc + l) ^^-^^ 

and, based on (5.5) we find ([24]) a similar, but independent result 

three of an infinite number of such recursive relationships ([10]). The latter two may be identified 
with known results by reversing the sums, giving 

.(.^. i-'r (^^ (2fc-i)c(2fc) A 

^^^"^ 2i-2" + 2n-2 (-7r2)'=r(2n-2A; + 2)y ' ^ ^ ^ 



thought to be new. 



^i^^i- 21-2" -2 I Z.(_^2)fcr(2n-2fc + l)) ^^-^^ 



, fe=o 



equivalent to [35], Eq.4.4(ll), and 



1 _ 21-2" 

equivalent to [35], Eq.4.4(9). 



.^o„^- n-^T fy.^ (l-2^^-^)C(2fc)(-l)n 

C(2^)-1_9i-2n 1 (27r)2fcr(2n-2fc + 2)J ^^"^^ 
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8.2 Novel series representations 

From (5.13), (5.16) and (5.27) with iV > 2, s = 2 find 



C(2) = -4a;^ (^arctan(-^) - l)j + Li2(e-2'^) - 2wlog(-l + e^-^) + 3^^ + 2a; (8.7) 

after identifying tfie fiypergeometric function with s = 2 in (5.16) from standard sources (e.g. [21], 
Chapter 9)^^. The series (8.7) is absolutely convergent by Gauss' test for arbitrary values of ui 
and reduces to an identity in the limiting case w — ^ 0. Similarly for the case s = 3 we find 



C(3) 



, m=l 



2 2 
TT m 



log(l + 



, ,2 



(8.8) 



,)- 1) +2a;Li2(e-2") + Li3(e-2-) 



-2w2iog(-l + e2") 
for n = 4 the result is 



(10a; + 3)w2 



^ m—l 



C(4) = ?'^' ( V(^arctan(— )-^ + J)|+2a;2Li2(e-2-)+2c.Li3(e-2-) 



TT^m^ 1 



+ Li4(e-2") - -Lj^logi-l + e^"^) 



and for n = 5 



(8.9) 



c(5)=^-M5:(. 



TT^ 1 
2 



■iog(i+ zi:::^) + + ^'^'Li2(e-^-) 



4 

r 



(8.10) 



+ 2a;2Li3(e-2-) + 2a;Li4(e-2-) + Li5(e-2-)-^a;4iog(-l + e2-) + ^ + ^. 

o 

When applied to (5.27), the general case follows from the well-known identifications ([21], 
Eq.(9.8.5)) 

TT arctan( ) 



2Fi(o> 1; oi -■ 



l/m^ = 



nm 



and 



2Fi(l, 1; 2; 



,)/m2 = ^ln(l + 



(8.11) 



B.12) 



By adding and subtracting the first n — 1 terms in the series representation of the hypergeometric 
function in (5.16), we find, for the case s = 2n + 2 



2Fi(l, n+ -; n+ -, 



( 



(n+o)(- 



^arctan(^)-i,g^^'"^' 

777, (J TT m. m ^ ^ — ^ 1 / . 



V 



7rm' 1/2 + k 

k=o ' 



'1 



^Compare the case N=0 with ([18], Eq. (42.1.1)) 



(8.13) 
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and for the case s = 2n + 3 



2 



Fi(l, n + 1; n + 2; 5-^) j = 



(8.14) 



1^ — , V _9 9 / 



^ln(l + ^^ ^> — ^^-^ 

k=0 



\ I 

for use in (5.13), (5.23) and (5.27). I can find no direct record in the literature, with particular 
reference to works involving Barnes double Gamma function (e.g. [7] ), that such representations 
are known, although, as will be shown below, they are equivalent to known results. 

8.3 Connection with multiple Gamma function 

From Section(5.2) with w = 1 we have, after some simplification, 

= 5 - 4 g <-"";;„7;<P">) + L.*-) - 21og,e^ - 1,. (8.15) 
Comparison of the above infinite series with [34], Eq.3.4(464) using z = i/-K identifies 



log ^ . , ; = 1— (6 + 41og(2 7r) + 2Li2(e-^) - 41og(e^ - 1) - 7r73) = 0.232662175652953 i 

\G(1 — j/Tr)/ 4 TT 

(8.16) 

where G is Barnes double Gamma function. From [34], Section 1.3, unnumbered equation 
following 26, we also find 

- mm = ^ - ^ (I (-1 

thereby recovering (8.7) after representing arctan in (8.7) as a logarithm. In the case s = 3, from 
(5.22) and (5.27) we have the series representation 

CO , . . , ^ (8-18) 
C(3) = 2Li2(e-2) + Li3(e-2) - 21og(e2 _ 1) + ^ _ 2 ^ ^ '^'"^ + 1" " "' 

m=l 

from which similar representions for the logarithm of the triple Gamma function can be 
obtained by recourse to [34], Eqs, 3.4(445) and (501). 

9 The Critical Strip 

In (4.17) and (4.18) with < (7 < 1 let 

f(s) = - ''"''^Y"^ (9.1) 
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Then, after adding and subtracting, we find a composition of symmetric and anti-symmetric 
integrals about the critical line cr = 1/2 : 

C(5) = /°°h(i) cos(pln(i)) sinh((l/2 - a) ln(t)) dt 

9{s) - f{s) Jo 

+ rm sin(pln(i)) cosh((l/2 - a) \n{t)) dt 

f{s) - g{s} Jo 

and simultaneously 



(9.4) 



C{s) = cos(pln(t)) cosh((l/2 - a) ln(t)) dt 

2if{s)g{s) 



(9.5) 



/ 

Jo 



h{t) sin(pln(t)) sinh((l/2 - a) \n{t)) dt . 



Let 



and 



m+9{s) 

pOO 

C{p)= h{t)cos{p\n{t))dt (9.6) 
Jo 

POO 

S(p)= / h{t)sm{pln{t))dt . (9.7) 
Jo 

For ({s) = on the critical line a = 1/2, both (9.4) and (9.5) can be simultaneously satisfied 
when both C{p) ^ and S{p) ^ only if 

1/(1/2 + zp)|2 = |.g(l/2 + zp)|2 (9.8) 

a true relationship that can be verified by direct computation. Set a = 1/2 in (9.4) and (9.5) and 
find 

C(l/2 + ^P)- + ip)+ 5(1/2 + i p) ^^^^ - /(1/2 + ip)-g{l/2 + i p) ^^'^ ^'"^^ 

Eq. (9.9) can be used to obtain a relationship between the real and imaginary components of 
C(l/2-|-v). Let 

C(l/2 + ip) = Cfl(l/2 + ip) + i C/(l/2 + ip). (9.10) 
Then, by breaking f {1/2 + ip) and 5(1/2 + i p) into their real and imaginary parts, (9.9) reduces to 

Cfi(l/2 + ip) _ C„(p)cos(pln(2 7r)) - C7p(p) sin(pln(2 7r)) 



C/(l/2 + ip) Cp{p) cos(pln(2 7r)) + C™(p) sin(pln(2 7r)) - 



(9.11) 



where 



Cp(p) =rK(l/2 + i/))cosh(7rp/2)+r/(l/2 + ip)sinh(7rp/2) (9.12) 
<^m(p) =rj(l/2 + ip)cosh(7rp/2)-rH(l/2 + ip)sinh(7rp/2) (9.13) 

in terms of the real and imaginary parts of r(l/2 + zp). This is a known result ([26]) that 

depends only on the functional equation (2.13)^^. In [26] it was shown that the numerator and 
denominator on the right-hand side of (9.11) can never vanish simultaneously, so this result 
defines the half-zeros^* (and usually brackets the full-zeros) of ({1/2 + I p) according to whether 



^^See the Appendix for further discussion on this point, 
^''only one of Cfl(l/2 -|- ip) or iJj(l/2 + ip) vanishes 



21 



the numerator or denominator vanishes. This was studied in [26]. The fuU-zeros^^ of ({1/2 + i p) 
occur when the following condition 

S{po) = 0, (9.14) 
C{po) = (9.15) 

is simultaneously satisfied for some value(s) of po- In that case, the defining equation (9.11) 
becomes 

Cr{1/2 + ip)' ^ _ C'{l/2 + ip)j ^ C^(p)cos(pln(2 7r)) - Cy(p) sin(pln(2 7r)) 

C/(l/2 + ip)' C'(l/2 + ip)fl Cp(p)cos(pln(2 7r)) + (7™(p)sin(pln(2 7r))-^/^ ^" ' 

where the derivatives are taken with respect to the variable p, and (9.16) is evaluated at p = po- 
Many solutions to (9.14) and (9.15) are known ([28]), so there is no doubt that these two 
functions can simultaneously vanish^^ at certain values of p. 

A further interesting representation can be obtained inside the critical strip by setting c = cr in 
(2.9) and separating ry(s) into its real and imaginary components 



rjR{a + pI) = ^a(^--^ 

(1 +t^)-'"/^ eparctan(t) (cos(0(t)) sin(7r a) coshj-K a t) -sm{n{t)) cosing) sinh(7r g f)) 
3 cosh^ {nat) — cos^ (tt a) 

(9.17) 

r?j(c7 + p7) = -la(i-) 

(1 + i2)-'^/2 e^'ar^Wt) (sin(0(i)) sin(7r cr) cosh(7r cri) + cos(0(t)) cos(7r cr) sinh(7rc7t))^^ 

3 cosh^ (tt a t) — cos^ (tt a) 

(9.18) 

where 

n{t) = -plog{a'^ (1 + 1'^)) + crarctan(t) . (9.19) 

In analogy to (9.4) and (9.5), the representations (9.17) and (9.18) include a term sin(7r(T) that 
vanishes on the critical line a = 1/2, both reducing to the equivalent of (3.7) in that casc^''. The 
integrands in both representations vanish as exp{—Tr at) as f — oo and are finite at t = 0; most of 
the numerical contribution to the integrals comes from the region < t = p. Numerical evaluation 
in the range <t can be acclerated by writing 

gparctan(t) _ (arctan(t)-7r/2) ^ (9.20) 

and extracting a factor e'' ^^"^ outside the integral in that range. It remains to be shown how the 
remaining integral conspires to cancel this factor in the asymptotic range p ^ oo. Finally, these 
two representations are invariant under the transformations p —p and t —t, demonstrating 
that r]{s) is self-conjugate, at least in the critical strip. 

10 Comments 

A number of related results that are not easily categorized, and which either do not appear in the 
reference literature or which appear incorrectly should be noted. 



2BCfl(l/2 + ipo) = C/(l/2 + ipo) = 

^^A proof that this cannot occur in (9.4) and (9.5) when a ^ 1/2 is left as an exercise for the reader. 
^'^It is suggested that such forms contain information about the distribution of zeros of r]{s) and hence C{s) inside 
the critical strip. 
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Comparison of (4.8) and (4.13) suggests the identification v = i. Such sums and integrals have 
been studied ([2], [7]) for the case s = —n in terms of Barnes double Gamma function G{z) and 
Clausen's function defined by 

ci.(.)=f:^, (10.1) 

fe=i 

with the (tacit) restriction that z & Di. That is not the case here. Thus straightforward 
application of that analysis will sometimes lead to incorrect results. For example [7], Eq.(4.3) 
gives (with n = 1) 

TT ^ cot(7rv)du = t log(2 tt) + log (10-2) 
further ([2], Eq.(13)) identifying 

To allow for the possibility that t G this should instead read 

fG(l-t)\ _^^_fM^t)\ , in^ , e-W%, 



valid for all t in terms of Bernoulli polynomials Bn{t) and polylogarithm functions. However, 
similar results are given in the equivalent, general form in [36], Eq. (5.12). Interestingly, both the 
Maple and Mathematica computer codes manage to evaluate^^ the integral in (10.2) in terms 
similar to (10.4) for a variety of values of n provided that the user restricts \t\ < 1, although the 
right-hand side provides the analytic continuation of the left when that restriction is lifted. This 
has application in the evaluation of C(3) according to the method of Section (8.3), where, using 
the methods of [34], Chapter 3, we identify 

x:ii^^=--^y^.^cot(..)d., (10.5) 

a generalization of integrals studied in [9]. Prom the computer evaluation of the integral, (10.5) 
eventually yields 

(10.6) 

y- C(2fc)t^'' ^l.^ Li,(e^-*) _ Li3(e^-VC(3) 

k=l 

as a new addition to the lists found in [34], Section (3.4) and [35], Sections (3.3) and (3.4), and a 
generalization of (5.22) when s = 3. 



A Appendix 

In [26], an attempt was made to show that the functional equation (2.13) and a minimal number 

of analytic properties of C(s) might suffice to show that C(,s) 7^ when a 7^ 1/2. This was 
motivated by a comment attributed to H. Hamburger ([19], page 50) to the effect that the 
functional equation for C,{s) perhaps characterizes it completely^^. The purported demonstration 

contradicts a result of Gauthier and Zeron ([16]), who show that infinitely many functions exist, 
arbitrarily close to C,{s)i that satisfy (2.13) and the same analytic properties that arc invoked in 

■^^Ncithcr [9], nor references contained therein, both of which study such integrals, Hst the computer results, which 
are likely based on [36], Eq. (5.12). 

■^^An unequivocal equivalent claim can also be found in [13]. As demonstrated by the counterexample given here, 
and in [16], this claim is false. 
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C = . (A.2) 



[26], but which vanish when ct ^ 1/2. ^° 

In fact, an example of such a function is easily found. Let 

v[s) = C sin(7r (s — sq)) sin(7r (s + sq)) sin(7r (s — Sq) sin(7r (s + Sg)) (A-l) 
where Sq is the complex conjugate of sq, So = co + is arbitrary complex, and C, given by 

1 

(cosh^(7r po) — cos^(7r ctq))^ 

is a constant chosen such that u{l) = 1. Then 

C.(s)^K5)C(s) (A.3) 

satisfies analyticity conditions (1-3) of [16] (which include the functional equation (2.13)). 
Additionally, Ci^{s) has zeros at the zeros of ({s) as well as complex zeros at s = ± (n + sq) and 
s = ± (n + Sq) and a pole at s = 1 with residue equal to unity, where n = 0, ±1, ±2, 

Thus the last few paragraphs of Section 2 of [26] beginning with the words "There are two cases 
where (2.9) fails... " cannot be vahd. However, the remainder of that paper, particularly those 
sections dealing with the properties of ^(1/2 + ip) as referenced here in Section 9 are valid and 
correct. I am grateful to Prof. Gauthier for discussions on this topic. 
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Abstract 

Abstract: Contour integral representations for Riemann's Zeta function and Dirichelet's 
Eta (alternating Zeta) function are presented and investigated. These representations flow 
naturally from methods developed in the 1800's, but somehow they do not appear in the 
standard reference summaries, textbooks or literature. Using these representations as a basis, 
alternate derivations of known series and integral representations for the Zeta and Eta 
function are obtained on a unified basis that differs from the textbook approach, and results 
are developed that appear to be new. 

1 Introduction 

Riemann's Zeta function C,{s), and its sibling, Dirichelet's (alternating zeta) function •q{s) play an 
important role in physics, complex analysis and number theory, and have been studied extensively 
for several centuries. In the same vein, the general importance of a contour integral representation 
for any function has been knovifn for almost two centuries, so it is surprising that contour integral 
representations for both C,{s) and ry(s) exist that cannot be found in any of the modern handbooks 
(NIST, [29], Section 25.5(iii); Abramowitz and Stegun, [1], Chapter 23), textbooks (Apostol, [5], 
Chapter 12; Olver, [30], Chapter 8.2; Titchmarsh, [37], Chapter 4; Whittaker and Watson, [39], 
Section (13.13)), summaries (Edwards, [12], Ivic, [19], Chapter 4; Patterson, [31]; Srivastava and 
Choi, [34] and [35]), compendia (Bateman, [14], Section 1.12 and Chapter 17), tables (Gradshteyn 
and Ryzhik, [17], 9.512; Prudnikov et. al, [32], Appendix II.7) and websites ([13], [29], [40], [41], 
[42], [43]) that summarize what is known about these functions. One can only conclude that such 
representations, generally discussed in the literature of the late 1800's, have been long-buried and 
their significance overlooked by modern scholars. 



It is the purpose of this work to disinter these representations, revisit and explore some of the 
consequences. In particular, it is possible to obtain series and integral representations of C(s) and 
77(5) that unify and generalize well-known results in various limits and combinations and 
reproduce results that are only now being discovered and investigated by alternate means. 
Additionally, it is possible to explore the properties of C,{s) in the complex plane and on the 
critical line, and obtain results that appear to be new. Many of the results obtained are disparate 
and difficult to categorize in a unified manner, but share the common theme that they are all 
somehow obtained from a study of the revived integral representations. That is the unifying 
theme of this work. To maintain a semblance of brevity, most of the derivations are only sketched, 
with citations sufficient enough to allow the reader to reproduce any result for her/himself. 

* mike@geometrics-unlimited.com 
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2 Contour integral representations for ({s) and t]{s) 

Throughout, I shall use s = a + ip, where s G € (complex), a,pG^R (reals) and n and A'' are 
positive integers. 

The Riemann Zeta function ({s) and the Dirichelet alternating zeta function ry(s) are well-known 
and defined by (convergent) series representations: 



C(s) = ^fc-^ a>l (2.1) 
fe=i 

oo 

r?(s)^-^(-l)'=fc-« a>0 (2.2) 



fe=i 



with 



,7(s) = (l-2i-)C(5). (2.3) 

It is easily found by an elementary application of the residue theorem, that the following 
reproduces (2.1) 

C(s) = 77 residue I t=fc (p —— — - dt a > 1 (2.4) 
^ Mk sin(7ri) 

and the following reproduces (2.2) 



oo „ 

r]{s) = n'^residue\t=k f cot(7rt) e'''*dt ct > (2.5) 

fe=i 

where the contour of each of the integrals encloses the pole at t = k. 

Convert the result (2.4) into a contour integral enclosing the positive integers on the t-axis in a 
clockwise direction, and, provided that u > 1 so that contributions from infinity vanish, the 
contour may then be opened such that it stretches vertically in the complex t-plane, giving 

C(s) = - / coi{TTt)dt 0<c<l,ceai (2.6) 

^ J c—ioo 

1 r°° 

= --/ {c + it)~^ coi{-K{c + it))dt a > 1, (2.7) 

2 J-oo 

a result that cannot be found in any of the modern reference works cited (in particular [36]). It 
should be noted that this procedure is hardly novel - it is often employed in an educational 
context (e.g. in the special case s = 2) to demonstrate the utility of complex integration to 
evaluate special sums ([34] Section 4.1), ([15], Lemma 2.1), and forms the basis for an entire 
branch of physics ([27]). 

A similar application to (2.5), with some trivial trigonometric simplification, gives 



. .c+ioo ^ 



r](s) = - / — — - dt (2.8) 

' 2y,_,^ sin(7rf) ^ ' 

+ dt. (2.9) 



1 /-°° (cH 

2 sin(7r 



{c + it)) 

Eq. (2.9) reduces, with c = 1/2, to a Jensen result ([20]), presented in 1895, and to the best of my 

knowledge, reproduced, again, but only in the special case c = 1/2, only once in the modern 
reference literature ([41]). In particular, although Lindelof's 1905 work [22] arrives at a result 
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equivalent to that obtained here, attributing the technique to Cauchy (1827), it is performed in a 
general context, so that the specific results (2.7) and (2.9) never appear, particularly in that 
section of the work devoted to C{s). In modern notation, for a meromorphic function f{z), 
Lindelof writes ([22], Eq. /// (4)) in general 

oo If 

E f{'^) = -7r-f^^°^i^ ^) /(^) dz, (2.10) 



7y— — oo 



where the contour of integration encloses all the singularities of the integrand, but never 
specifically identifies f{z) = z~^. Much later, after the integral has been split into two by the 
invocation of identities for cot{'K z), he docs make this identification ([22], Eq. (6)), arriving 

(with 2.3) at (3.7 - see below) a now-well-known result, thereby bypassing the contour integral 
representation (2.7). In a similar vein, Olver ([30], page 290) reproduces (2.10), employing it to 
evaluate remainder terms, the Abcl-Plana formula and a Jensen result, but again, never writes 
(2.7) explicitly. Perhaps the omission of an explicit statement of (2.7) and (2.9) from Chapter IV 
of Lindelof ([22]) is why these forms seem to have vanished from the historical record, although 
citations to Jensen [20] and Lindelof [22] abound in modern summaries^. 

At this point, it is worthwhile to quote several related integral representations that define some of 
the relevant and important functions of complex analysis, starting with Riemann's original 1859 
contour integral representation of the zeta function ([33]) 

where the contour of integration encloses the negative t-axis, looping from t = — oo — zO to 
t = —oo + iO enclosing the point t=0. 

The equivalence of (2.11) and (2.1) is well-described in texts (e.g. [35]) and is obtained by 
reducing three components of the contour in (2.11). A different analysis is possible however (e.g. 
[12], Section 1.6) - open and translate the contour in (2.11) such that it lies vertically to the right 
of the origin in the complex t-plane^ with a <0 and evaluate the residues of each of the poles of 
the integrand lying a.t t = ± 2nni to find 

oo 

C(s) = — r(l - s)(27r«)^(exp(i7rs/2) - exp(-i7rs/2)) V k"''^. (2.12) 
This can be further reduced, giving the well-known functional equation for ({s) 

= (2,3) 

valid for all s by analytic continuation from a < 0. This demonstrates that (2.11) and (2.6) are 
fundamentally different representations of C{s), since further deformations of (2.11) do not lead to 
(2.7). 

Ephemerally, it has been claimed [41] that the following complex integral representation (missing 
a factor 27r in [41]) 



C(^) = T77r^/ dt a>0 (2.14) 



2{s-l)J_^ cosh2(7ri) 

can be obtained by the invocation of the Cauchy-Schlomilch transformation [4]. In fact, the 
representation (2.14) follows by the simple expedient of integrating (2.7) by parts in the special 
case c = 1/2 (see (3.16) below), as also noted in [41]. Another result 



^Recently, Srivastava and Choi [35], pages 169-172 have reproduced Lindelof's analysis, also failing to arrive 
explicitly at (2.7); although they do write a general form of (2.8) ([35], Section 2.3, Eq. (38)), it is studied only for 
the case c = 1/2. 

^and thereby vanishes 
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(i_2i-«)r(s + i)yo cosh2(t) 

obtained in [4] on the basis of the Cauchy-Schlomilch transformation does not follow from (2.14) 

as erroneously claimed in [41]. 

Based on (2.11), a number of other representations of ({s) are well-known. From the tables cited, 
the following are worth noting here: 
([14], Eq.l.l2(6)) 



([14], Eq.l.l2(4)) 

r(s) Jo exp(t) - 1 



1 roo 

C(s) = — - / — — -dt C7 > 1 , (2.17) 



C(s) = T^FTT^-TT / —TYT^'^t. a>l (2.18) 



and ([17], Eq. 3.527(1)) 

1 (2a)(''+i) f°° t 

4 r(s + 1) Jo ^M?(ai) 

By way of comparison, Laplace's representation ([17], Eq. 8.315(2)) 

,(c+it) 



1 _ 1 r 

T{s)~2^iJ_, 



:dt (2.19) 



{c + ity 

defines the inverse Gamma function. 

3 Simple reductions 
3.1 The case < c < 1 

The simplest reductions of (2.7) and (2.9) are obtained by uniting the two halves of their range 
corresponding to f < and < > 0. With the use of 

{-{c + i ty + {c-i ty) i = 2 (c2 + f2)(V2) sin(s arctan(t/c)) (3.1) 

and 

{c + ity + {c-i ty = 2 (c^ + i2)(«/2) cos(s arctan(i/c)) (3.2) 
one finds (with < c < 1 and < c < 1 respectively) 

(3.3) 

cos(s arctan(t/c)) cos(7r c) sin(7r c) — sinh(7r t) cosh(7r t) sin(s arctan(f/c)) 



(is) = -[ 
Jo 



(c2 +^2)^/2 (cosh2(7ri) - cos2(7r c)) 

a > 1 



-dt 



and 

^ .o.^/o^ ^2. o. ' ' ^-^dt, (3.4) 



-sin(sarctan(t/c)) cos(7rc) sinh(7rt) + cos(s arctan(f/c)) sin(7rc) cosh(7r<) 
(c2 + t^yli (cosh2(7rt) - cos2(7rc)) 



the latter valid for all s. In the case that c = |, (3.3) reduces to 

sinh(./)^u(.Hrc-|.a.(2f)) 
' Jo (l + 4i2)V2cosh(7rt) ' ^ ^ 
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equivalent to the well-known result"^ (which also appears as an exercise in many textbooks - e.g. 
[35], Eq.2.3(41)) 

C(a)-^-2^ /°° sin(^arctan(20)e-* 
^^'> - s-1 ^ io (1 + 4 t^y/' cosh(7r t) "'^ • 

Similarly, (3.4) reduces to the known ([22], IV.III(6)) result 

io (l + i2)V2cosh(7rt/2)"'*- ^"^-^^ 
Comparison of (3.5) and (3.6) immediately yields the known result ([17], Eq.3. 975(2)) 

^ 2(^-1) /"^ sin(sarctan(2i))e'^* , , 

^(^) = T37 + 2Vo (1 + 4.^)^2^4^) '^- ^'-'^ 

Notice that (3.6) and (3.7) are valid for all s. Simple substitution of c ^ in (3.4) gives 

??(s) = -sin(s7r/2) / — - — -dt a < (3.9) 



/o sinh(7r t) 
which, after taking (2.3) into account, leads to 



sin(s TT /2) 

(21-" - 1) Jq sinh(7ri) 



CG^) = i^^/ -TT^dt, a<0 (3.10) 



displaying the existence of the trivial zeros of ({s) at s = —2N, since the integral in (3.10) is 
clearly finite at each of these values of s. The result (3.10) is equivalent to the known result (2.16) 
after taking (2.13) into account. Other possibilities abound. For example, set c = 1/3 in (3.3) to 
obtain 

_ os-i (2 siii(s arctan(.)) sinh(2 tt t /3) - V3 cos(.s arctan(.))) 
' Jo (l+t2)«/2(2cosh(2^</3) + l) ■ ^' ' 

Further, a novel representation can be found by noting that (2.9) is invariant under the 
transformation of variables t ^ t + p together with the choice c = a along with the requirement 
that < o" < 1. This effectively allows one to choose; c = ,s, yielding the following complex 
representation for ri{s), as well as the the corresponding representation for ({s) because of (2.3), 
valid in the critical strip < u < 1: 



' r <• + "'"• (3.13) 



2 y_ooSin(7r (s + it)) 

Because (3.1) and (3.2) are valid for c G £, it follows that (3.4) with c = s is also a valid 
representation for r]{s), as well as the the corresponding representation for ({s) with recourse to 
(2.3), again only valid in the critical strip < < 1. Another interesting variation is the choice 
c = 1/p in (3.4), valid for p > 1, again applicable to the corresponding representation for ({s) 
because of (2.3). Such a representation may find application in the analysis of ^(0" + ip) in the 
asymptotic limit p ^ 00. Finally, the choice c = a, 0<a<l leads to an interesting variation 
that is discussed in Section 9. 



^attributed to Jensen by Lindelof, and a special case of Hermite's result for the Hurwitz zeta function ([39], 
Section 13.2). 
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3.2 Simple Recursion via Partial Differentiation 

As a function of c, ({s) and r]{s) are constant, as are their integral representations over a range 
corresponding to < N — c < 1; however the integrals change discontinuously at each new value 
of c that corresponds to a new value of TV. Thus, the operator ^ acting on (2.7) and (2.9) 
vanishes, except at non-negative integers c = N, when it becomes indefinite. Therefore, it is of 
interest to investigate this operation in the neighbourhood of c = N/2, as well as the limit N. 
Applied to (2.9)^, the requirement that 

dr]{s)/dc = (3.14) 

and identification of one of the resulting terms, yields 

2s J _^ sm (7r(c + tt)) 
a new result with improved convergence at infinity. Similarly^ , operating with 4- on (2.7) yields 



a.+i) = ^r .y/t' (3.16) 

2s y_ooSm^(7r {c + it}) 

Using (3.1) and (3.2) to reduce the integration range of (3.15), gives a fairly lengthy result for 
general values of c. If c = 1/2, (3.15) becomes 

2-. ^-(l + .^)-V^sinh(.t/2)sin(.arctan(t)) 
' s Jo cosh2(7ri/2) 
which, by letting s ^ s — 1 can be rewritten as 

,(. + l) = ^ii^M£)+?!:iZ[ rcos(sarctan(t))(l + .^)-V^sinh(.V2)^^^^ 
s s Jo cosh (7rt/2) 

and, with regard to (2.3) 



2(s-l)(1^2-^ 
^ s(l-2-«) 

2^-1 TT f°° cos(s arctan(t)) (1 + i2)-V2 sinh(7r ^2) 



(3.19) 



7 

Jo 



tdt . 



s {1-2-') Jo cosh^(7rt/2) 
Similarly^ , in the case that c = 1/2, (3.16) reduces to the well-known result ([22], Eq. 4,//7(5)'^) 

Cis + 1) = ^ r ^os{s^rctan(f^ 

' s Jo il + f^y/^cosh\nt/2)'^^- ^^-^"^ 

If c in (3.15)^ , one finds, a finite result^. 



, S TT 



7rsin( — ) .oo^l-. w ^^ 

V{s)= ^ ; / —r-j-^dt, a<0 3.21 

s-1 Jo smh (tt*) 



^and equivalent to integration by parts 

■'It should be noted that Srivastava, Adamchik and Glasscr ([36]) study an incomplete form of (3.16) in order to 
obtain representations for sums both involving C('^)> s-nd CC") by itself, witliout commenting on the significance of 
the contour integral with infinite bounds (3.16). 

®(3.18) and (3.19) are the first two entries in a family of recursive relationships that will be discussed elsewhere 
(in preparation). 

'^attributed to Jensen [20] 

^For the case c = in (3.16) see the following section 
^also equivalent to integrating (3.9) by parts 
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and, when c = 1, one finds the unusual forms 

iirr r{l + rie' ^)(-*) cosh(7r re'^)e'^ 



r TT f 

s J I sinh''(7rri) 



(3.22) 



sinh^(7rr e'^) 
r TT (1 + f2)(-s/2) (,og(-g arctan(r t)) cosh(7rrt) 



d6» 



or, 



r]{s + 1) = 1 + —- / ^ ^ '- de 

r TT /■°° (1 + t2)(-s/2) (,Qg(g arctan(r <)) cosh(7r r t) 



(3.23) 



sinh (-Trri) 



■dt 



valid for any^" < r < 1 and all s, after deformation of the contour of integration in (3.15) either 
to the left or right of the singularity at t = 1, respectively. 

4 Special cases 
4.1 The case c=0 

With respect to (3.16), setting c = together with some simplification gives 



C(S) = -TT^ 



,sin(7rs/2) t^-' , „ 



s-1 Jo sinh^(t) 

reproducing (2.18), with a= 1, after taking (2.13) into account. This same case (c = 0) naturally 
suggests that the integral in (2.6) be broken into two parts 

C(s) = Co(s) + Ci(s) (4.2) 

where 

Co(s) = ^ (f t~'coi(nt)dt , (4.3) 



2 



(^i(s) = - / cot{-Kt) dt + - cot{-Kt) dt 



(4.4) 



such that the contour in (4.3) bypasses the origin to the right. It is easy to show that the sum of 
the two integrals in (4.4) reduces to 

/oo 
coth.{-Kt)dt a>l (4.5) 

sin(7rs/2) . , f°° e^-""'^ , 

= — — + sin(7r s/2) / —— — - dt , (4.6) 
s - 1 ^ ' ' Ji sinh(7r t) ' ^ ' 

the latter valid for all s. One, of several commonly used possibilities of reducing (4.3) is to apply 
([18] (54.3.2)) 

2 °° 

t cot(7rt) = ^t^'' C(2fc) \t\ < 1 (4.7) 

giving 



^"particularly r = l/7rorr = l/2 
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fc=0 

convergent for all s. Such sums appear frequently in the literature (e.g. [8], [34] Section (3.2)), 
and have been studied extensively when s = —n ([34], [7]). In the case s = 2n, (4.8) reduces to 
Co(s) = C{s) and (4.6) vanishes. 

Along with (4.1) the integral in (3.16) can be similarly split, using the convergent series 
representation 

^ oo 

= (2^-1) k) v' \v\ < 1 (4.9) 

sin (ttv) ^ 

where the Bernoulli numbers B2k appearing in the original ([18], Eq.(50.6.10)) are related to 
C(2fc) by 



Thus the representation 



1 in r t^-" , 7rsin7r.s/2) t^-* , 

C(s) = — -/ — ^ dt i '—^ 5 dt (4.11) 

' 2 s-1 J-i sin2 TTt s-1 Ji sinh'(7rt) ^ ' 



becomes (formally) 



^ s-1 [t^ fr'o 2A;-s A sinh2(7ri) y ^ ' 

after straightforward integration of (4.9) with a < 0. In (4.12), the integral converges, and the 
formal sum conditionally converges or diverges according to several tests. Howver, further analysis 
is still possible. Following the method of [34], Section (3.3) applied to (4.8), together with [34], 
Eq. 3.4(15), we have the following general result for \v\ < 1 

fe=o -"J fc=l 

1 1 r 

= — + -v'' t-" {^{l+t)-^il-t))dt a<2 (4.14) 

2 s 2 Jq 

so that (formally), after premultiplying (4.14) with 1/v, then operating with v^-^, 



E C(2fe)(2 fc 1)^^^ =-^ + h' f^'-' (*'( 1 + *) + 1 - t)) dt (4.15) 



2k-s 2s 2 

k=0 



= -^ + lv' ft^-' (tt' (1 + coe{nt)) -hdt. a<2 (4.16) 

Notice that the right-hand sides of (4.14) provides an analytic continuation of the left-hand side 
for \v\ > 1 and the right-hand side of (4.16) gives a regularization of the left-hand side. Putting 

together (4.12) and (4.16), incorporating the first term of (4.16) into the second term of (4.12), 
setting V = i and constraining cr > eventually yields a new integral representation^^ 

-k' 'sm(— ) .oo ^ ^ 



icompaxe with [38], Eq. (2.8). 
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valid for < (7 < 2, a non-overlapping analytic continuation of (4.1) that spans the critical strip. 
Additionally, using (2.13) followed by replacing s ^ 1 — s gives a second representation 



valid for — 1 < ct < 1, a non-overlapping analytic continuation of (2.18). Over the critical strip 
< (7 < 1 representations (4.17) and (4.18) are simultaneously true and augment the rather short 
list of other representations that share this property ([6], as cited in [14]^^). See Section 9 for 
further application of this result. 

4.2 The case c > 1 

Interesting results can be obtained by translating the contour in (2.6) and (2.9) and adding the 
residues so-omitted. From (2.1) and (2.2) define the partial sums 

N 



Civ(s)^^A;-« (4.19) 
fe=i 

N 

w(s) = -l](-l)'^fc-^ (4.20) 



fe=i 

let cn,N = 1, 2, 3, ... with 

N<CN<N + 1 (4.21) 

and adjust^^ all the above results accordingly by redefining c—^Cn. This gives, for example, the 
generalization of (2.6) 

C(s) = Cjv(s) - i^Sn,c^ - - / t-' cotiwt) dt ,a>l (4.22) 

^ JcN-ioo 

and (2.9) 

nis) = V.is) + i^^.,.. + - y^^_^^ dt (4.23) 

so (3.3) and (3.4) with cn become integral representations of the remainder of the partial 
sums C(s) — C,n{s) and r}{s) — r]N{s), the former corresponding to a special case of the Hurwitz 
zeta function ({s, N). 

In the case that cn = N = 1, one finds 



and 



^, , 1 /■°°cosh(7ri)sin(sarctan(i)) , ^ 

_ s + l ^J-e--^^^^^^ ^^^^^ 



2(s-l) Jo (l-Fi2)V2sinh(7ri) 



, , 1 /'°° sin(sarctan(t)) ,^ ,^ 

with the latter two results valid for all s. Neither of (4.24) nor (4.26) appear in references; (4.25) 
corresponds to a limiting case of Hermite's representation of the Hurwitz zcta hmction ([14], Eq. 



l^unnumbered equations immediately following 1.12(1), also reproduced in [35], Eqs. 2.3(44)-(46). 

^^If cjv = ^ an extra half-residue is removed and the contour deformed to pass to the right of the point t = N. 
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1.10(7)). In analogy to the developments leading to (3.8), equations (4.24), (4.26) and (2.3), 
together with the following modified form of ([17], Eq.(9.51.5)) 

C(,)- '^"'^ I /■- sin(.arctan(^))(l + e-*) 

(2.-l)(s-l) + 2«-iyo (l + i2)V2sinh(7rt) ^ ' 

can be used to obtain the possibly new result 

s-3 /"^ sin(sarctan(i))e''* , , 

For c = 1/2 + corresponding to general partial sums, two possibilities are apparent. Based on 
(3.1), (3.2), (3.3) and (3.4)1^ we find 



^11 r 



sinh(7r i (1 + 2 A^)) (1 + t2)-s/2 sin(s arctan(f)) , 

at (T > 1 

cosh2(7rt (1/2 + 7V)) 

(4.29) 

(l/2 + 7V)i-« /■'^(l + i2)-V2sin(sarctan(t)) 



s-1 

and 



which together extend (3.5) and (3.7) respectively. Application of (3.15) and (3.16) with (3.1) and 
(3.2), also with c = 1/2 + A'', gives the remainders 

as + -J1I1±^ r (l + ^-);^-cos(^arctan(0) 

fe^ ^ -^0 cosh2(7rt(l/2 + 7V)) ^ ^ 



k=l 

and 



''(^ + i) + Eyr 
fc=i 



(4.33) 



TT (-1)^ (1/2 + Ny-' f°° sinh(7r t (1/2 + TV)) (1 + i2)-«/2 sin(s arctan(t)) 



/ 

Jo 



cosh2(7rt(l/2 + AT)) 

valid for all s. For the case c = N applying (4.22) and (4.23) gives the remainders 

C(.) - f 1 = - J- + N^- rcosh(.Art)(l + .^)-V^sin(.arctan(t)) 

2Ar* Jo sinh(7rAri) ^ ^ 

= -^ + ^+2Ar-^ r (l+^^)-/^sin(.arctan(^)) 

2iV- + s-l+ io e2^-*-l ^ 

and 

^ A;^ 2Ar« ^ ^ Jo sinh(7rAri) ^ ^ 

which extend (4.24) and (4.26), the latter valid for all s. It should be noted that (4.35), valid for 
s G € generalizes well-established results (e.g. [17], Eq. 3.411(12)) that are only valid in the 



*for comparison, see also [30], Section 8.3, and [12], Section (6.4) and Chapter 7. 
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limiting case s = n^^. Omitting the integral term from (4.30) gives an upper bound for the 
remainder; similarly, omitting the integral term in (4.35) gives a lower bound for the remainder, 
provided that ,s G in both cases. In addition, many of the above, in particular (4.30) and (4.35) 
appear to be new and should be compared to classical results such as [19], Theorem 1.8, page 23 
and pages 99-100. 



5 Series Representations 
5.1 By Parts 

By simple manipulation of the integral representations given, it is possible to obtain new series 
representations. Starting from the integral representation (3.10), integrate by parts, giving 

({s)= sm(7rs/2)7r^ 



(2(i--)-l)(-s/2+l) 
^^-iF.(-| + l; ^,2-|; ^)cosh(7rt) 



dt cr < 



sinh^(7r t) 

where the hypergeometric function arises because of [23], Eq. 6. 2. 11(6). Write the hypergeometric 
function as a (uniformly convergent) series, thereby permitting the sum and integral to be 
interchanged, apply a second integration by parts, and find 



( \ °° 2k foc ,l-s+2k 

JT^) FpIT2)i l^-^' 

Identify the integral in (5.2) using (2.18) giving 

, _ 2^ w'-^ sin(7r s/2) ^ 2-^>' ({-s + 1 + 2 k)r{-s + 2 + 2k) 
(,{s) - A. r(2fc + 2) ^^""^^ 

valid for all values of a by the ratio test and the principle of analytic continuation. This result, 
which appears to be new, can be rewritten by extracting the fc = term, and applying (2.13) to 
yield 

, _ 2^7r^-isin(7r5/2) ^ 2-^1^ Cjs + 1 + 2 k)Ti-s + 2 + 2 k) 

21--2 + S ^ r(2fc + 2) ' ^ ' 

fc=i ^ ' 

a representation, valid for all values of s, that is reminiscent of the following similar sum recently 
obtained by Tyagi and Holm ([38] Eq.(3.5)) 



7r°-- (2-2«-2'=)C(-s-|-l + 2fc)r(-s-|-l + 2fc) 
Cis) = T3^sm(. s/2) ^ i r(2fc + 2) ' ^>'' 

(5.5) 



r(2fc + 2) 



Eq. (5.4) can also be rewritten in the form 



I 16 



1 f. 2-^^C(^ + 2fc)r(l + 5 + 2fc) 

^~r(s)(2--l-s)^ r(2A; + 2) ^^'^^ 



-'^Although [17], Eq. 3.411(21) gives (with a missing minus sign) an integral representation for finite sums such 
as those appearing in (4.35), specified to be valid only for integer exponents (fe"), that particular result in [17] also 
appears to apply after the generalization n — > s. 

^^by Gauss' test 
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by replacing s — 1 — s and applying (2.13). 

With reference to (2.7), express the numerator cosine factor as a difference of exponentials each of 
which in turn is written as a convergent power series in tt (c + it), recognize that the resulting 
series can be interchanged with the integral and each term identified using (2.9) to find 

^ /^V^ (-^V4)'^??(.s--2fc)\ 

A more interesting variation of (5.7) can be obtained by applying (2.3) and explicitly identifying 
the fc = term of the sum giving 

a convergent representation if ct > 1 (by Gauss' test). An immediate consequence is the 
identification 

C(2) = -y C(0). (5.9) 

because only the k = 1 term in the sum contributes when s = 2. Additionally, when s = 2n, (5.8) 

terminates at n terms, becoming a recursion for C(2't.) in terms of C(27^ — 2fc), k = 1, . . .n - see 
section 8 for further discussion. If the functional equation (2.13) is used to extend its region of 
convergence, (5.8) becomes 

.y^ C(l-s + 2fc) (2^-^^-1-1) 
r(2A: + l)r(s-2fc) 

^(-)=2 "\2-^-l)cos(../2) « 

which augments similar series appearing in [34], Sections 4.2 and 4.3, and corresponds to a 
combination of entries in [18] (Eqs. (54.6.3) and (54.6.4)). 

5.2 By splitting 
5.2.1 The lower range 

After applying (2.13), the integral in (4.1) is frequently (e.g. [35]) and conveniently split in two, 
defining 

p UJ J.S 

C-{s)= -TjjrJt (5-11) 

^0 smh (t) 

and 



with 



Tis + l) 
2 

for arbitrary w. Similar to the derivation of (4.12), after application of (2.13) we have 



C(s)4l^=C-W + C+W (5.13) 



k=0 

which can be rewritten using (4.10) as 



^""h - 1/ 2) r(fc + 1) (2 fc - 1 + s) 



(5.14) 
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~ {-'^) r(t + i)c(2fc) 



U) < TV 



(5.15) 



Substitute (2.1) in (5.15), interchange the order of summation, and identify the resulting series to 
obtain 

= ^i^H^^ ^ 2 + 2^2 + 2^ -;^)-^-th(a;) + ^ , (5.16) 

a convergent representation that is vahd for all s. 

Other representations are easily obtained by applying any of the standard hypergeometric linear 
transforms as in (6.4) - see below. One such leads to 



^2 2^^r(fc+i + |)^^V2mVa;2 + l 



f+'- w^coth(w) + 



SUJ 



s-l 



1 



(5.17) 



a rapidly converging representation for small values of w, valid for all s. The inner sum of (5.17) 
can be written in several ways. By expanding the denominator as a series in l/{t^rin?) with 
r < 2p — 1 and interchanging the two series we find the general form 



E 

m=0 



T{p + m) (-l)™f-2m-2p^(-2^ + 2j9-r) 
r(p) r(m + 1) 



(5.18) 



where the right hand side continues the left if \t\ > 1 and r — 2p> —1. Alternatively, [18], Eq. 
(6.1.32) identifies 



E 

m=0 



1 1 TT -uf^y\ 

(m2 t^ + y^) " V ^ 27^^^°^ 



(5.19) 



so for p = fc + 1 and r = in (5.18), each term of the inner series (5.17) can be obtained 
analytically by taking the fc*'' derivative of the right-hand side with respect to the variable at 
y = 1, leading to a rational polynomial in coth(a;) and sinh(a;). The first few terms are given 
below for two choices of w. 



The first few values of «;(0, A; + 1, tt^/w^) in (5.18) 



analytic 



0J = 1 



u = n 



-1/2+ wcoth(a;)/2 



1^1 cv^ 
2 4 sinh^w 



H — w coth(a;) 
4 



3a; j.u3 1 iv2 , cjcothcjf — 3+2li; ) 

coth w — i coth uj H , y ^ 

I 8+3 ui^ 

~^ If) sinh^c^' 



0.156518 1.076674 
0.927424 X 10"^ 0.306837 

0.795491 X 10"'^ 0.1.34296 



Application of (5.18) to (5.17) leads to a less rapidly converging representation that can be 
expressed in terms of more familiar functions 
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. , ,^ (5-20) 



^2 2'^^^^ r(fc+| + |)r(m + l) 

— w^cothfw) -I . 

s — 1 

A second transform of the intermediate hypergeometric equation yields another useful result 
2sL0^-' ^ ^=+1)^ u.'^ „^(_l)™^(m + A;+i + |)C(2m + 2A; + 2) 

^-^'^-r(i/2 + 5/2)^/^^ ^/^^ r(m + i)r(fc + i)(2fc + i + s) 

„ , ,s-l 

— a;*coth(a;) 



s- 1 

(5.21) 

Both (5.20) and (5.21 ) can be rc-ordcrcd along a diagonal of the double sum-'^^, allowing one of the 
resulting (terminating) series to be evaluated in closed form. In either case, both eventually give 



oo 



2 



C-(«) = -2sa;-i^(-J^) ^(^^-^^coth(a;), (5.22) 

reducing to (4.8) in the case cj = tt taking (2.13) into consideration. Furthermore, by applying 
(2.1) and interchanging the resulting series, a further (more slowly converging) representation is 
found 

s — 1 ^ 1 -1 2 

m=l 

which differs from (5.16)^^. Series of the form (5.22) have been studied extensively ([34], [35]) 
when s = n for special values of w. The representations given above will find application in 
Section 8. 



5.2.2 The upper range 

Rewrite (5.12) in exponential form and expand the factor (1 — e~^*)~^ in a convergent power 
series, interchange the sum and integral and recognize that the resulting integral is a generalized 
exponential integral ^s{z) usually defined in terms of incomplete Gamma functions ([25], [29], 
Section 8.19 ) as follows 

Y,,{z) ^ z-'-^T {\- s,z) (5.24) 
giving a rapidly converging series representation 

oo 

C+(s) =4w"+^ ^/cE_^(2wfc). (5.25) 

k=l 

A useful variant of (5.25) is obtained by manipulating the well-known recurrence relation ([29], 
Eq. 8.19.12) for E^(0) to obtain 

E_,(z) = ^ (e-^(s + z) + s{s - 1) E2_«(z)) . (5.26) 

Apply (5.26) repeatedly A^" times to (5.25) and find 

"Er=o E^=o k) = Er=o E^=o k - m) 

^*and could have also been obtained by recourse to contiguity relations. 
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c,(.) = iog(e- - 1) + 2..- ^ + r(. + i)g 

2cj-''-^-^r(.s + l) ^ Ejv+2-.(2fca;) 
r(s - TV - 1) 2^+1^ fc^+i 

(5.27) 

in terms of logarithms and polylogarithms {Lik(z)), a result that diverges numerically as A'' 
increases, but which also, usefully, terminates if s = n, N > n > 0. This is discussed in more 
detail in Section 8. 

Another useful representation emerges from (5.25) by writing the fimction Es(..) in terms of 
confluent hypergeometric functions. From [25], Eq.(2.6b) and [21], Eq.(9.10.3) we find 

(5-28) 

C+(s) = 4a;«+i^fcE_,(2A;w) = 

fc=i 

A S-\- 1 ^ 

2i-« r(s + 1) C(s) - V fc 6-2'=- iFi(l; 2 + s;2ku), 

^ + ^ fe=i 

where the series on the right hand side converges absolutely if cr > 1 independent of u>. 
Comparison of (5.25), (5.13) and (5.28) identifies 

4 S+ 1 ^ 

C-(s) = — V&e-2'='"iFi(l; 2 + s; 2koj) (5.29) 

s + 1 

k=l 

which can be interpreted as a transformation of sums, or an analytic continuation, for the various 
guises of Q-{s) given in Section(5.2.1), if both, or either, of the two sides converge for some choice 
of the variables s and/or cj, rcspcictively. The fact that a term containing C,{s) naturally appears 
in the expression (5.28) for C+(s) suggests that there will be a severe cancellation of digits 
between the two terms C-(s) and (^+{s) if (5.13) is used to attempt a numerical evaluation of C,{s), 
as has been reported elsewhere ([11]). Essentially, (5.13) is numerically useless (with common 
choices of arithmetic precision) for any choices of p > 10, even if w is carefully chosen in an 
attempt to balance the cancellation of digits. However, with 20 digits of precision, using 20 terms 
in the series, it is possible to find the first zero of C(l/2 + ip) correct to 10 digits, deteriorating 
rapidly thereafter as p increases. For any choice of oj , the convergence properties of the two sums 
representing C,-{s) and C+(s) anti-correlate, as can be seen from the Table in the previous Section. 
This can also be seen by rewriting (5.28) to yield the combined representation 

2i-r(s + l)C(s)=4u;«+i (^A:(E_,(2fcu;) + i^^iii±li^^±^i^^)j . (5.30) 

In (5.30), the first inner term (E_s(2fca;)) vanishes like exp(— 2 fca;)/(2 fccj) whereas the second 
term (ii^i/(s + 1)) vanishes like l/fc(^+*) — exp(— 2 /c a;)/(2 few) so for large values of k cancellation 
of digits is bound to occur 

It should be noted that (5.28) is representative of a family of related transformations. Although 
setting s = in (5.28) does not lead to a new result, operating first with ^ followed by setting 
s = leads to a transformation of a related sum 

(5.31) 



a;arctanh(^^) ^a;arctanh(— |^) i ^ i ^ 
^Ei(2fcw) = r7=4^+9log(-) + ?77 + ^ 



k=l ^ k=2 



^^In fact, the coefficients of both asymptotic series prefaced by the exponential terms cancel exa<;tly term by term. 
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where the first term corresponding to fc = 1 has been isolated, since it contains the singularity 

_d_ 

ds 



that occurs due to divergence of the left-hand sum if lj = iw. Similarly, operating with 4- on 



(5.28) at s = 1, leads to the identification 

oo 

E + ^) = -1 - l°g(2) + log(e2 - 1) (5.32) 
fc=i ^ 

(which can alternatively be obtained by writing the infinite sum as a product). A second 
application of ^ on (5.28), again at s = 1 leads to the following transformation of sums 



fe=i 



(5.33) 



1^, ,7r2p + 1 , ^ 1 V-2 + 2e% ^Ei(2A;) 

k=l k=l 

This transformation can be further reduced by applying an integral representation given in [29] , 
Section 6.2, leading to 

Correcting a misprinted expression given in ([3], Eq. (4.8), [34] and [35], Eq. 3.3(43)) which 
should read 

^ ^ C(2 k) = log(7r csc(7r V^))—, (5.35) 
followed by a change of variables, eventually yields the representation 

f: ^'""-,'r-""' = -^/i(>OE(ie' - ie-') - M*)).*. (5.36) 

Substituting (5.32), (5.34)) and (5.36) into (5.33) followed by some simplification, finally gives 

.OS ( '^-'-y-"l ^ = + ^ + .(1) - i .»S-^(2) + 2 . (5.37) 



/ 

Jo 



2v J V 12 2 ' ' 2 



a result that appears to be new. 

5.3 Another Interpretation 

These expressions for the split representation can be interpreted differently. For the moment, let 
a < —3, and insert the contour integral representation [25], Eq. (2.6a) for E_s(z) into the 
converging series representation (5.25), noting that contributions from infinity vanish, giving 

a(.) = -4c.-«f::^ r"°"'-"'""''d. (5.38) 

s+i + t 

where the contour is defined by the real parameter c < —2. With this caveat, the scries in (5.38) 
converges, allowing the integration and summation to be interchanged, leading to the identification 

2 TT I Jc-i oo S+l+t 

The contour may now be deformed to enclose the singularity at t = — s — 1, translated to the right 
to pick up residues of r(— f) at t = n, n = 0, 1 . . . and the residue of C(~l — t) at t = —2. After 
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evaluation, each of the residues so obtained cancels one of the terms in the series expression (5.15) 
for (s) leaving (after a change of integration variables) 

r(. + i)C(.)_ 1 / r(. + i)CW 



the standard Cauchy representation of the left-hand side, valid for all s, where the contour 
encloses the simple pole at t = s. 



6 Integration by parts 

Integration by parts applied to previous results yields other representations. 
Prom (2.17), we obtain^" 

1 f°° 

'^('^ = - Y{s-l) I i'"'log(l-exp(-i))rft a>l. (6.1) 



and from (4.24) we find 

1 s /•°° log(sinh(7rt)) cos((s + 1) arctan(^)) 
^^'^=2-nJo (l + ^^)a+^)/^ ^^-'^ 

as well as the following intriguing result^^ 



1 1 V^^"^(^)r(^-1) /■o° tsiu(.sairUm(^)) 
C(sj =0 + TsZ + ^ / ■ .2/ 21*1(1/2, s/2; 3/2; ) dt 

2 2 p(_) 7o smh {-Kt) 



(6.3) 



/■°°fcosh(7rt)cos(sarctan(t)) ^ ,^ . „. „, , 

-W ■ u ^.2^ 2F1 1/2, a/2; 3/2; -t^ . a>0 

Jo smh(7rt) (1 +t^) 

This result has the interesting property that it is the only integral representation of which I am 
aware in which the independent variable s does not appear as an exponent. Many variations of 
(6.3) can be obtained through the use of any of the transformations of the hypergeometric 
function, for example 

2Fi(l/2, s/2; 3/2; -t^) = (1 + t')-^/^ 2Fi(l, s/2; 3/2; f/il + f)), (6.4) 
which choice happens to restore the variable s to its usual role as an exponent. 



7 Conversion to Integral Representations 

Ftom (5.4), identify the product C(...)r(...) in the form of an integral representation from one of 
the primary definitions of ({s) ([17], Eq. (9.511)) and interchange the scries and sum. The sums 
can now be explicitly evaluated, leaving an integral representation valid over an important range 
of s 



. / 2ssinh(V2) , , , A 

4^ sin(. s/ 2) r ' V 1 - ' + ) 

2 - 21+" + 5 2" Jo e* - 1 ^ ' 

Each of the terms in (7.1) correspond to known integrals ([17], Eqs. (3.552.1) and (8.312.2)), and 
if the identification is taken to completion, (7.1) reduces to the identity ({s) = ({s) with the help 



^"known to Maple, but only in the form of a Mellin transform 
discovered by the Maple computer code 
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of (2.13). Other useful representations are immediately available by evaluating the integrals in 
(7.1) corresponding to pairs of terms in the large parenthesis, giving 



7r^2«-2 /^sin(7rs)2^ /■°°* '(-l + ^o^M^)) l\ 
^(-)= (l-2^)cos(.a/2) |-^yo e^-1 ''^W)) ^'''^ 

and 

({s) = tt'-'^ sin(7r s/2) ( f u""-^ (w/sinh(t;) - 1) e"" dv - — — ) . a<2 . 

^ ^ ^ ' \Jo \ ! I sm(7rs)r(l + s); 

(7.3) 

8 The case s = n 

Several of the results obtained in previous sections reduce to helpful representations when s = n. 

8.1 Finite series representations 

From (5.4) we find, for even values of the argument, 

" (-l)fc7r^fc(2n-2fc^l)C(2n-2fc) 
ih r(2/c + 2) 

'^('-) = -2- F^^TTTT^l ' 

from (5.8) we find 

1 ^ (-7r^)^(l-2i-^"+^^)C(2n-2fc) 
C(2«)- 2(2-2" _ 1)1. r(2fc + l) ^^-^^ 

and, based on (5.5) we find ([24]) a similar, but independent result 

three of an infinite number of such recursive relationships ([10]). The latter two may be identified 
with known results by reversing the sums, giving 

.(.^. i-'r (^^ (2fc-i)c(2fc) A 

^^^"^ 2i-2" + 2n-2 (-7r2)'=r(2n-2A; + 2)y ' ^ ^ ^ 



thought to be new. 



^i^^i- 21-2" -2 I Z.(_^2)fcr(2n-2fc + l)) ^^-^^ 



, fe=o 



equivalent to [35], Eq.4.4(ll), and 



1 _ 21-2" 

equivalent to [35], Eq.4.4(9). 



.^o„^- n-^T fy.^ (l-2^^-^)C(2fc)(-l)n 

C(2^)-1_9i-2n 1 (27r)2fcr(2n-2fc + 2)J ^^"^^ 
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8.2 Novel series representations 

From (5.13), (5.16) and (5.27) with iV > 2, s = 2 find 



C(2) = -4a;^ (^arctan(-^) - l)j + Li2(e-2'^) - 2wlog(-l + e^-^) + 3^^ + 2a; (8.7) 

after identifying tfie fiypergeometric function with s = 2 in (5.16) from standard sources (e.g. [21], 
Chapter 9)^^. The series (8.7) is absolutely convergent by Gauss' test for arbitrary values of ui 
and reduces to an identity in the limiting case w — ^ 0. Similarly for the case s = 3 we find 



C(3) 



, m=l 



2 2 
TT m 



log(l + 



, ,2 



(8.8) 



,)- 1) +2a;Li2(e-2") + Li3(e-2-) 



-2w2iog(-l + e2") 
for n = 4 the result is 



(10a; + 3)w2 



^ m—l 



C(4) = ?'^' ( V(^arctan(— )-^ + J)|+2a;2Li2(e-2-)+2c.Li3(e-2-) 



TT^m^ 1 



+ Li4(e-2") - -Lj^logi-l + e^"^) 



and for n = 5 



(8.9) 



c(5)=^-M5:(. 



TT^ 1 
2 



■iog(i+ zi:::^) + + ^'^'Li2(e-^-) 



4 

r 



(8.10) 



+ 2a;2Li3(e-2-) + 2a;Li4(e-2-) + Li5(e-2-)-^a;4iog(-l + e2-) + ^ + ^. 

o 

When applied to (5.27), the general case follows from the well-known identifications ([21], 
Eq.(9.8.5)) 

TT arctan( ) 



2Fi(o> 1; oi -■ 



l/m^ = 



nm 



and 



2Fi(l, 1; 2; 



,)/m2 = ^ln(l + 



(8.11) 



B.12) 



By adding and subtracting the first n — 1 terms in the series representation of the hypergeometric 
function in (5.16), we find, for the case s = 2n + 2 



2Fi(l, n+ -; n+ -, 



( 



(n+o)(- 



^arctan(^)-i,g^^'"^' 

777, (J TT m. m ^ ^ — ^ 1 / . 



V 



7rm' 1/2 + k 

k=o ' 



'1 



^Compare the case N=0 with ([18], Eq. (42.1.1)) 



(8.13) 
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and for the case s = 2n + 3 



2 



Fi(l, n + 1; n + 2; 5-^) j = 



(8.14) 



1^ — , V _9 9 / 



^ln(l + ^^ ^> — ^^-^ 

k=0 



\ I 

for use in (5.13), (5.23) and (5.27). I can find no direct record in the literature, with particular 
reference to works involving Barnes double Gamma function (e.g. [7] ), that such representations 
are known, although, as will be shown below, they are equivalent to known results. 

8.3 Connection with multiple Gamma function 

From Section(5.2) with w = 1 we have, after some simplification, 

= 5 - 4 g <-"";;„7;<P">) + L.*-) - 21og,e^ - 1,. (8.15) 
Comparison of the above infinite series with [34], Eq.3.4(464) using z = i/-K identifies 



log ^ . , ; = 1— (6 + 41og(2 7r) + 2Li2(e-^) - 41og(e^ - 1) - 7r73) = 0.232662175652953 i 

\G(1 — j/Tr)/ 4 TT 

(8.16) 

where G is Barnes double Gamma function. From [34], Section 1.3, unnumbered equation 
following 26, we also find 

- mm = ^ - ^ (I (-1 

thereby recovering (8.7) after representing arctan in (8.7) as a logarithm. In the case s = 3, from 
(5.22) and (5.27) we have the series representation 

CO , . . , ^ (8-18) 
C(3) = 2Li2(e-2) + Li3(e-2) - 21og(e2 _ 1) + ^ _ 2 ^ ^ '^'"^ + 1" " "' 

m=l 

from which similar representions for the logarithm of the triple Gamma function can be 
obtained by recourse to [34], Eqs, 3.4(445) and (501). 

9 The Critical Strip 

In (4.17) and (4.18) with < (7 < 1 let 

f(s) = - ''"''^Y"^ (9.1) 
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Then, after adding and subtracting, we find a composition of symmetric and anti-symmetric 
integrals about the critical line cr = 1/2 : 

C(5) = /°°h(i) cos(pln(i)) sinh((l/2 - a) ln(t)) dt 

9{s) - f{s) Jo 

+ rm sin(pln(i)) cosh((l/2 - a) \n{t)) dt 

f{s) - g{s} Jo 

and simultaneously 



(9.4) 



C{s) = cos(pln(t)) cosh((l/2 - a) ln(t)) dt 

2if{s)g{s) 



(9.5) 



/ 

Jo 



h{t) sin(pln(t)) sinh((l/2 - a) \n{t)) dt . 



Let 



and 



m+9{s) 

pOO 

C{p)= h{t)cos{p\n{t))dt (9.6) 
Jo 

POO 

S(p)= / h{t)sm{pln{t))dt . (9.7) 
Jo 

For ({s) = on the critical line a = 1/2, both (9.4) and (9.5) can be simultaneously satisfied 
when both C{p) ^ and S{p) ^ only if 

1/(1/2 + zp)|2 = |.g(l/2 + zp)|2 (9.8) 

a true relationship that can be verified by direct computation. Set a = 1/2 in (9.4) and (9.5) and 
find 

C(l/2 + ^P)- + ip)+ 5(1/2 + i p) ^^^^ - /(1/2 + ip)-g{l/2 + i p) ^^'^ ^'"^^ 

Eq. (9.9) can be used to obtain a relationship between the real and imaginary components of 
C(l/2-|-v). Let 

C(l/2 + ip) = Cfl(l/2 + ip) + i C/(l/2 + ip). (9.10) 
Then, by breaking f {1/2 + ip) and 5(1/2 + i p) into their real and imaginary parts, (9.9) reduces to 

Cfi(l/2 + ip) _ C„(p)cos(pln(2 7r)) - C7p(p) sin(pln(2 7r)) 



C/(l/2 + ip) Cp{p) cos(pln(2 7r)) + C™(p) sin(pln(2 7r)) - 



(9.11) 



where 



Cp(p) =rK(l/2 + i/))cosh(7rp/2)+r/(l/2 + ip)sinh(7rp/2) (9.12) 
<^m(p) =rj(l/2 + ip)cosh(7rp/2)-rH(l/2 + ip)sinh(7rp/2) (9.13) 

in terms of the real and imaginary parts of r(l/2 + zp). This is a known result ([26]) that 

depends only on the functional equation (2.13)^^. In [26] it was shown that the numerator and 
denominator on the right-hand side of (9.11) can never vanish simultaneously, so this result 
defines the half-zeros^* (and usually brackets the full-zeros) of ({1/2 + I p) according to whether 



^^See the Appendix for further discussion on this point, 
^''only one of Cfl(l/2 -|- ip) or iJj(l/2 + ip) vanishes 
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the numerator or denominator vanishes. This was studied in [26]. The fuU-zeros^^ of ({1/2 + i p) 
occur when the following condition 

S{po) = 0, (9.14) 
C{po) = (9.15) 

is simultaneously satisfied for some value(s) of po- In that case, the defining equation (9.11) 
becomes 

Cr{1/2 + ip)' ^ _ C'{l/2 + ip)j ^ C^(p)cos(pln(2 7r)) - Cy(p) sin(pln(2 7r)) 

C/(l/2 + ip)' C'(l/2 + ip)fl Cp(p)cos(pln(2 7r)) + (7™(p)sin(pln(2 7r))-^/^ ^" ' 

where the derivatives are taken with respect to the variable p, and (9.16) is evaluated at p = po- 
Many solutions to (9.14) and (9.15) are known ([28]), so there is no doubt that these two 
functions can simultaneously vanish^^ at certain values of p. 

A further interesting representation can be obtained inside the critical strip by setting c = cr in 
(2.9) and separating ry(s) into its real and imaginary components 



rjR{a + pI) = ^a(^--^ 

(1 +t^)-'"/^ eparctan(t) (cos(0(t)) sin(7r a) coshj-K a t) -sm{n{t)) cosing) sinh(7r g f)) 
3 cosh^ {nat) — cos^ (tt a) 

(9.17) 

r?j(c7 + p7) = -la(i-) 

(1 + i2)-'^/2 e^'ar^Wt) (sin(0(i)) sin(7r cr) cosh(7r cri) + cos(0(t)) cos(7r cr) sinh(7rc7t))^^ 

3 cosh^ (tt a t) — cos^ (tt a) 

(9.18) 

where 

n{t) = -plog{a'^ (1 + 1'^)) + crarctan(t) . (9.19) 

In analogy to (9.4) and (9.5), the representations (9.17) and (9.18) include a term sin(7r(T) that 
vanishes on the critical line a = 1/2, both reducing to the equivalent of (3.7) in that casc^''. The 
integrands in both representations vanish as exp{—Tr at) as f — oo and are finite at t = 0; most of 
the numerical contribution to the integrals comes from the region < t = p. Numerical evaluation 
in the range <t can be acclerated by writing 

gparctan(t) _ (arctan(t)-7r/2) ^ (9.20) 

and extracting a factor e'' ^^"^ outside the integral in that range. It remains to be shown how the 
remaining integral conspires to cancel this factor in the asymptotic range p ^ oo. Finally, these 
two representations are invariant under the transformations p —p and t —t, demonstrating 
that r]{s) is self-conjugate, at least in the critical strip. 

10 Comments 

A number of related results that are not easily categorized, and which either do not appear in the 
reference literature or which appear incorrectly should be noted. 



2BCfl(l/2 + ipo) = C/(l/2 + ipo) = 

^^A proof that this cannot occur in (9.4) and (9.5) when a ^ 1/2 is left as an exercise for the reader. 
^'^It is suggested that such forms contain information about the distribution of zeros of r]{s) and hence C{s) inside 
the critical strip. 
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Comparison of (4.8) and (4.13) suggests the identification v = i. Such sums and integrals have 
been studied ([2], [7]) for the case s = —n in terms of Barnes double Gamma function G{z) and 
Clausen's function defined by 

ci.(.)=f:^, (10.1) 

fe=i 

with the (tacit) restriction that z & Di. That is not the case here. Thus straightforward 
application of that analysis will sometimes lead to incorrect results. For example [7], Eq.(4.3) 
gives (with n = 1) 

TT ^ cot(7rv)du = t log(2 tt) + log (10-2) 
further ([2], Eq.(13)) identifying 

To allow for the possibility that t G this should instead read 

fG(l-t)\ _^^_fM^t)\ , in^ , e-W%, 



valid for all t in terms of Bernoulli polynomials Bn{t) and polylogarithm functions. However, 
similar results are given in the equivalent, general form in [36], Eq. (5.12). Interestingly, both the 
Maple and Mathematica computer codes manage to evaluate^^ the integral in (10.2) in terms 
similar to (10.4) for a variety of values of n provided that the user restricts \t\ < 1, although the 
right-hand side provides the analytic continuation of the left when that restriction is lifted. This 
has application in the evaluation of C(3) according to the method of Section (8.3), where, using 
the methods of [34], Chapter 3, we identify 

x:ii^^=--^y^.^cot(..)d., (10.5) 

a generalization of integrals studied in [9]. Prom the computer evaluation of the integral, (10.5) 
eventually yields 

(10.6) 

y- C(2fc)t^'' ^l.^ Li,(e^-*) _ Li3(e^-VC(3) 

k=l 

as a new addition to the lists found in [34], Section (3.4) and [35], Sections (3.3) and (3.4), and a 
generalization of (5.22) when s = 3. 



A Appendix 

In [26], an attempt was made to show that the functional equation (2.13) and a minimal number 

of analytic properties of C(s) might suffice to show that C(,s) 7^ when a 7^ 1/2. This was 
motivated by a comment attributed to H. Hamburger ([19], page 50) to the effect that the 
functional equation for C,{s) perhaps characterizes it completely^^. The purported demonstration 

contradicts a result of Gauthier and Zeron ([16]), who show that infinitely many functions exist, 
arbitrarily close to C,{s)i that satisfy (2.13) and the same analytic properties that arc invoked in 

■^^Ncithcr [9], nor references contained therein, both of which study such integrals, Hst the computer results, which 
are likely based on [36], Eq. (5.12). 

■^^An unequivocal equivalent claim can also be found in [13]. As demonstrated by the counterexample given here, 
and in [16], this claim is false. 
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C = . (A.2) 



[26], but which vanish when ct ^ 1/2. ^° 

In fact, an example of such a function is easily found. Let 

v[s) = C sin(7r (s — sq)) sin(7r (s + sq)) sin(7r (s — Sq) sin(7r (s + Sg)) (A-l) 
where Sq is the complex conjugate of sq, So = co + is arbitrary complex, and C, given by 

1 

(cosh^(7r po) — cos^(7r ctq))^ 

is a constant chosen such that u{l) = 1. Then 

C.(s)^K5)C(s) (A.3) 

satisfies analyticity conditions (1-3) of [16] (which include the functional equation (2.13)). 
Additionally, Ci^{s) has zeros at the zeros of ({s) as well as complex zeros at s = ± (n + sq) and 
s = ± (n + Sq) and a pole at s = 1 with residue equal to unity, where n = 0, ±1, ±2, 

Thus the last few paragraphs of Section 2 of [26] beginning with the words "There are two cases 
where (2.9) fails... " cannot be vahd. However, the remainder of that paper, particularly those 
sections dealing with the properties of ^(1/2 + ip) as referenced here in Section 9 are valid and 
correct. I am grateful to Prof. Gauthier for discussions on this topic. 
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